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AEST8ACT 



This st'ady rsvisws the history and development of serine 
electric propulsion drives, the types of electric propul- 
sion, and the inherent losses which occur within the synch- 
ronous AC machines typically used for high-power propulsion 
systems. 

A thorough review of the literature pertaining to heat 
transfer in electrical machinery is made. In particular, 
the use of liquid cooling in various flow configurations, 
including buoyan cy- driven t her mosyphons and two-phase ther- 
losyphcns, is analyzed. 

Forced- liqu id cooling is feasible, but the required 
rotating seals are a problem in reliability. Closed-loop 
thermosyphcn cooling appears feasible at high rotational 
speeds, athough a secondary heat exchange through the shaft 
is required. Closed, two-phase thermosyphons and heat pipes 
are also feasible, but require forced'-air circulation for 
heat rejection to the ambient. Since all of these concepts 
deserve additional attention, areas for further research and 
development are recommended. 
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SOHENCIATORE 



Latin Synbcls 

a tube radius (m) 

Ac K fi2/g 

Cp specific heat (kJ/kg.K) 
d tube diamerer (m) 
f friction factor 

condensation factor for equation (3.6) 

F stress ratio = x, / 6g o. 
g acceler a'^ion of gravity (m/s2) 

Gr Grashof number = g g ^Td^/ 

Grr Rotational Grashcf number = H 2 g ^Td^/ 
hfg heat of vaporization (kJ/kg) 

H radius cf rotation (m) 
j vclumetric flux (m^/s) 

i* dimensicnless volumetric flux = j p ^ (go (p^ -p ) ) - 1/ 2 

y g y 

J Rctational Reynolds number = Sld2/ v 

k thermal conductivity (W/m.K) 

Ktcr-critical transition value (for equations 2.11-2.15) 

Kl = JRe (for equation 2. 14) 

Kt = J2/Re (for equations 2.11 & 2.12) 

L length cf tubs (m) 

p pressure (Pa) 

P pressure (Pa) 

Er Prandtl number 

g specific heat flux (W/m2) 

Q Heat load (W) 

Ra Rayleigh number (Gr*Pr) 

Far Rctational Rayleigh number (Grr*Pr) 

Fe Reynclds number = Vd/v 

1 0 



to Hcs^ty number = V / y)d 
S s'.iirl number = 1/ Ro 

t temperature (K or °C) 

V velocity (m/s) 

Greek Symbols 

a void fraction or thermal diffusiviry (unitless or la^/s) 
S ccefficient of thermal expansion (VK) 

5 film thickness (m) 

V ratic of specific heats 

A difference in values 

T temperature gradient (K/m) 

p density (kg/m^) 

u visccsity (N.s/m^) 

a surface tension (K/m) 

V kinematic viscosity (m^/s) 

0 . angular velocity (1/s) 

Subscripts 

0 stationary reference 

am arithmetic mean temperature difference 
b bulk ccnditicn 

E boiling section 

c ceded section 

cr critical value 

E entrainment value 

h heated section 

1 laminar conditions 

1 liquid state 

s scnic valu* (J1=1) 

sat saturation value 

t turbulent conditions 

vapor state 

1 1 



V 



w 



m 



wall condition 
maan condition 
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I. IN7RODOCTIOH 



The larine engineer and naval architect are charged with 
creating a design of a vessel from a broad set of require- 
ments for the ship system prepared by the customer. This 
requires a very structured procedure for preparing th® 
design-based coraparisions and economic trade-off srudies at 
each and every stage cf the development of the design. The 
end result is a complete set of design requirements from 
which to construct the final product (this result should 
meet the original design requirements!). 

This structured procedure has the form of the well-known 
"design spiral" [Ref. 1 ], which is shown in Figure 1.1 for 
the propulsion machinery. This procedure encompasses the 
total ship design, including (but not limited to): 

liSln F- CPU Is ion System 
Shaft horsepower 
Propeller RPM 

Specific fuel consumption and bunker capacity 
Space and weight objectives 
Adaptability to ship configuration 

A uxili arv Shin Sy stems 
Power and lighting 

Steam-galley, deck, and heating systems 
Heating, ventilating, and air conditioning 
Firefighting and ballasting 
Fresh water 
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Encin gq ring S vst am s 
k r. c ho r ha n d i ir. g 

St-csring engine and bridge teiemetering ccntrol 
Cargo handling gear 
Crane systems 

El ectr o nic a nd Ma viaat i on S ys tems 

Communication, exterior and interior 
Radar 

Loran, Omega, etc., navigation aids 

Military electronics, sensors, command and control 
systems, weapons directors, tactical data systems, 
and electronic countermeasures. 

In order to design the vessel within economic constraints, 
or to use the life cycle cost as the function to optimize, 
it is imperative to define initially the devices used to 
provide the power for the vessel. They must provide most for 
the least, in terms of the fuel that they expend and the 
volume and weight that they possess. 

The selection of the main propulsion plant requires 
matching the power of the generating device with a transmis- 
sion (in most cases) , the propulsor, other ship systems and 
the hull. If one limits the choices of propulsors to fixed- 
pitch propellers and controllable and reversible- pitch 
propellers, which are currently realistic for large vessels, 
the number of possible permutations remains quite sizable, 
as shown in Figure 1.2. The horsepower and RPM requirements 
of modern warships require the use of a transmission to 
couple the high RPM of an economical prime mover with the 
low RPM required at the propulsor for high propulsive 
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MISSION 




Figur 



e 




n Spiral 



efficiency. This transmission may be either mechanical or 
elactrical. The mechanical reduction gear has been 
exploited in the past to minimize volume/space requirements. 
When coupled to a reversible prime mover, such as a steam 
turbine, it has teen proven to be quite acceptable. 
However, the electrical transmission has several attractive 
features, such as ease of propulsor speed and direction 
control, flexibility in the the number and location of 
controls, freedom of installation and machinery layout, 
ccuplina of various prime movers independently to propulsor 
drives (as well as ether auxiliary applications) , and th® 
ease of auxiliary thruster installation for maneuvering 
control. 

The choice of the propulsion— plant includes many impor- 
tant factors, such as: 



1. Reliability 

2. Maintainability 

3. Space and arrangement requirements 
u. Weight requirements 

5. Type of fuel required 
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6. Fusl consuniction 

• Fractional power and ■transient perforinance 

8. Interrelations with auxiliaries 

9. Reversing capability 

10. Operating personnel 

11. Rating limitations 

12. Costs 
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Figure 1.2 Propulsion Plant ilternatives, 



These factors have gone through many cyclic changes 
based on the fluctuating economy. The recent emphasis on 
fuel costs and availability has driven the fuel 
consumption/prime mover efficiency to a much higher level of 
concern than it was previously given. The cost of acquisi- 
tion of new vessels has skyrocketed# as well as the costs of 
maintenance and repair. These escalating costs are a chal- 
lenge to the engineer of today. The economic considerations 
of the initial cost of a plant# the cost of maintenance and 
repair# and the cost to operate a plant from its 
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construction to its rerirstnent, also coined ”lif-? c^cla” 
cost, is the driving force that controls the design of the 
lodern warship. 

Of the list of prime movers shown in Figure 1.2, the 
lif 5-cycle cost of the marine gas turbine is one of the 
lowest. The economy of operation is maximized within a 
rather narrow range of RPM's. This has been tolerated in 

the current employment of this device, since the cost of 
utilizing an electric transmission to accomplish the goal of 
operating the prime mover within its maximum economy range 
did not offset the loss of efficiency experienced by the gas 
turbine. This was due to the technology that existed within 
the field of electrical controls at the rime. The advent of 
electronic frequency control of the output of the generator, 
independent of RPM or poles, and the frequency control of 
the synchronous motor, independent of its supply frequency, 
makes the use of the electrical transmission much more 
attractive. Although the initial cost will be higher, the 
COSO of operation, maintenance, and repair of the system 
will rapidly offset the initial cost differential [Ref. 2], 
[Ref. 3], and [Ref. 4]. An analysis of the economic and 
engineering considerations involved with this mechanical vs. 
electrical drive selection was made by Eric Gott and S. 0. 
Svensson [Ref. 5] in 1974; they clearly illustrated the 
nature of the problem and possible solutions. In their 
paper, it was evident that the selection of some of the 
possibly more appropriate forms of electric drives was 
discouraged due to the economics of the choice. High-power 
frequency-controlling equipment was very expensive and 
bulky. Recent advances in high-power, solid-state equipment 
have now again promoted these alternatives into viable 
contention. 
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A. aiSTOSICAL DEVELOPS ENT 

From th^= sarly attempts by the Russians in the 1890 »s -c 
power small passenger launches, and the British and rrnnch 
experiments with small subm ersibles , to the present, the use 
of electric propulsion has made a slow evolution. The O.S. 
Navy began with the Jupiter in 1911. it was a 55 00 hp/shaft 
vessel equipped with wound rotor induction motors, which 
fulfilled a 30~ysar lifespan (terminated by surface warfare 
activity in 1943, while serving as an aircraft carrier 
(Lar.gely)) . This led to several major programs, summarized 
helo w : 

* AC drive installed on five battle ships of 
30,000 + hp (1918-1 922) 

* AC drive applied to carriers Saratoga and Lexington of 
200,000hp (1920) 

* Synchronous ac drive introduced on four Coast Guard 
cutters (19 19) 

* Synchronous ac drive in the liner Normandie at 

160,000hp ( 1935) 

* Synchronous ac drive and full dynamic braking in 

cutters (19 40) 

In the O.S. Navy Bureau of Ships Manual (BDSHIPS 
250-660-2), dated 1 March 1945, the text states that the 
Navy's inventory of ships that used electric drive (in addi- 
tion to above) was "....about two hundred destroyer escorts 
and about one hundred other ships including tankers, cargo 
transports, troop transports, and store ships." All of these 
applications were with steam turbines (which were speed- 
adjustable and reversible) and eventually this design lest 
favor to the less expensive mechanical reduction gear 
systems which were somewhat lighter. The weight of these 
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(3600SPii) and twenty-eight to eighty pole moxors 
(256-9CBEM) . The large number of poles required in moxor the 
for speed reduction caused the motors to be extremely larce 
and heavy. The advent of the use of synchronous motors 
reduced the weight slightly since the motor experienced much 
fewer slip losses. A good historical background on the 
early stages of electric propulsion is found in the book 
written by Commander S.M. Robinson, USN [Ref. 6]. 

Current applications in the O.S. Navy are shown in 
Table I. 



TABLE I 





O.S. Navy 


Electric Drive 


Ships 


Type 


1 Shies 


SHP 


^ Motors 


AS 


2 


15,000 


1 


AS 


2 


11,520 


8 


TARC 


1 


10,000 


2* 


TAGOS 


1 


1,600 


2 ** 


♦ also has U- 


1200 hp thrusters 




** also has 1- 


550 hp thr 


ust er 





E. ELECTRIC PROPULSICH COHCEPTS 

An area of current research in the field of rotating 
electric equipment that shows considerable promise with 
regard to electric propiiision is the use of acyclic (homo- 
polar) superconductive DC devices. These devices produce 
exceptionally high torque and horsepower in a very small 
volume and weight. This is of great value in advanced 
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hullforai propulsion. i large aiaount of research in this 
area to cstsrmne ccnf igur ation, peak casualty potential, 
maintainability of the required cryogenic equioment., as well 
as the adaptibility of the systems for this cryogenic appa- 
ratus must be performed before actual implementation of this 
configuration can occur in a specific design. Thus, for 
rapid implementation in a short-term project (3 — 5yrs) , the 
choices remaining in electric propulsion to be considered 
are direct-curren t (DC) and alternating-current (AC) 

devices. The categcry of AC devices is further subdivided 
into inductive and synchronous devices. All three devices 
have teen used for marine propulsion and each has ixs cwn 
inherent advantages and disadvantages. 

1 . CC Mac h ine ry 

The DC electric drive system has been used exten- 
sively in the past, primarily in low to medium power ranges 
(1000 - 6000 hp/shaft) with extremes of 400 hp/shaft and 

19,600 hp/shaft noted. This employment of DC drives 
normally • utilizes 1-4 prime movers powering one or more 
motors per shaft (duplicating the system if more than one 
shaft exists) . 

Primary advantages of this system include; 

1. Ease of Cont rol . Control being effected by varying 
the generator voltage through field control-a process 
that lends itself well to remote location(s). 

2- gql tip le Control S tatio ns. The ease of which the 
generator voltage and polarity of the motor connections 
(either field or_exciter) are varied allows multiple 
control stations to be designed for flexible ship 
control in tight maneuvering situations. 

3. Ad aptab il ity to va rin g pro peller- hul l cha rac t eris - 
tics. The propeller-hull charact er ist ics change 

2 1 



d~2. !n £"t ic 2 . 1 iy w h €n t-hs v-^sssl is "towinij or ico~br£ 2 ,’*'inc. 
Tha ability to match th£ prime mov^r (rur.iiir.a at 
cptiiTium to the cropeiler thr'^iioh the ’D~or.’'’ is 

motor and still provide maximum torgus is a soecific 
advantage of this system, without excess engine or 
6l-=ctrical capacity, through a vide range of propell‘=r 
BFMs. 

Primary disadvantages include: 

1. Limited ^wer Range. The top-end range of approxi- 
mately 10,000 hp/shaft limits the vessel size that can 
utilize this system. 

2. Generator S peed L imi tations . The maximum speed 
that the generator can be operated, due to problems 
primarily associated with commutation, is typically 
limited to under 1000 BPM. This limits the selection 
of the prime mover if one doss not utilize mechanical 
reduction gearing to match the characteristics. Figure 
1.3 shows the normal relationship of RPM vs. KW. 

3- Height and S ize . The weight of the DC system 
increases dramatically with power rating as shewn in 
Figure 1.4 (this figure shows weight in pounds/SH?) ; 
note that total weight becomes almost linear with SHP 
at higher SHP values. This precludes high power appli- 
cations of conventional DC systems above 10,000 
hp/shaft ( acycl ic-homo polar devices are excluded from, 
this analysis) . The physical size of the device also 
has this dramatic functional relationship to rating as 
shown in Figure 1.5. 

From the above discussion of advantages and disad- 
vantages, the consideration of DC electric propulsion for 
use in a warship design that would require on the order of 
40,000 hp and utilize gas turbines in their most efficient 
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Figure 1 .3 



Variation of DC Generator Speed 
with Bating. 




Figure 1.4 



Variation of DC Generator Weight 
with Eating. 



operating range 
of this paper 



is net realistic. Therefore, the remainder 
will only consider the AC systems. The 
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SHAFT HOSSEPOWER 



Figure 1.5 Variation of DC Generator Size 

with Bating. 

inherent losses associated with these DC systems though are 
essentially the same as those in the AC systems. Figures 
1.1 through 1.5 are from '‘Marine Engineering", R. L. 
Harringtcn, editor. 

2 . ^ machi nerv 

The devices included in the conventional AC 
machinery area fall into twc categories: induction and 

synchronous. The principal difference between the twc 
systems is in the motor application. The induction motor 
system provides exceptional torque at low RPM and maximum 
slip { 200 %). This high torque response at low RPM makes the 
induction motor system preferable in some applications. For 
naval applications, this system is not as efficient as the 
synchronous system by a few critical percentage points. The 
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induction motor hi s losses that are a function of the slip 
incurred bs-.ween rctor and tha stator field. Ihis slip 

iTt2Ei3.i.rid.ncr Xoss^s 5 .^^' 



nha synchronous motor in 



is rslatsd to the RPS and lead. Th 
induction motor are the same as 
general. Another critical difference between rhe two 
systems is the power factor difference. In the induction 
motor, the rotor current is induced from the stator field. 
This induced field that is created as a result of this 
current is the mechanism from which the mechanical rotation 
is derived. Since this current is induced from the stator 
(supplied current), it causes the curr=>nt load seen by the 
supplying device to ’’lag" the supplying voltage and there- 
fore cause the power factor (a measure of the amount of lag 
or lead) *o be less than unity. A power factor of unity 
indicates that the load current and the load voltage are in 
phase and maximum real energy transfer is occurring. The 
induction motor is normally operated wish a lagging power 
factor of approximately 80/t [Eef. 6] and an overall effi- 
ciency cf approximately 90%. In the synchronous system, the 
power factor is normally unity and is controllable . by the 
external' motor rotor excitation which results in the synch- 
ronous system's power factor being adjustable and a unity 
power factor being easily obtained. Note here also that the 
power factor in this system can be made to actually "lead" 
and thereby correct a system "lag" to other loads that might 
exist in the fully integrated system. The advantages cf AC 
systems in general are: 

1. Hig h Ef fic iency. Losses in the synchronous system 
are approximately 6% (including excitation losses) as 
compared to the induction motor system losses that 
range from approximately 8 to 10%. 

2. Fle xib il it y of In stallation . The direct drive 
motor can b® installed in the vessel to minimize shaft 
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ler.gth to the propeller {sines physically smaller than 
DC motors) and the generator can be installed vinhout 
rsducrion gearing and operated direcrly coupled to rh'^ 
prime mover. 

3- Dual Ose of Fropul s ion Power. The AC system allows 
for the power from propulsion units to be used for 
c^her functions when not required for propulsion; this 
leads to the creation of a fully-integrated sysrem. 

ATailable in larqe P ow er Rating s . The conventional 
AC systems can be built in any power rating required 
with a reasonable upper limit of approximately 60,000 
hp. The larger system's size is being reduced by 
raising the system voltage (which is not limited by the 
commutation requirements of DC sysrems) . The system 
voltages typically in use are from 2,300 to about 7,500 
volts. These voltages are easily transformed for other 
uses in addition to propulsion. 

The selection of the main propulsion system there- 
. fore is not simple. The design environment must be specified 
and the prime mover is selected. Once these parameters are 
defined, the appropriate selection of transmission, control, 
and propulscr can be chosen. For the purpose of the propul- 
sion cf a vessel of approximately 8,000-ton displacement and 
the mission of a normal warship of the destroyer/cruiser 
variety using the marine gas turbine as the prime mover and 
life-cycle costs the design criteria, the selection of the 
AC synchronous drive (possibly fully integrated) appears as 
a very appropriate choice. Some of the design configura- 
tions possible within this choice, and the associated 
parfcrmance of each, are shown in Figure 1.6 as presented in 
[Ref. 7]. Additional possible configurations are shown in 
Figures 1.7 and 1.8 . The end result of the implementation 

c-^ these technologies would result in substantial savings 
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utilizes 'th'r 



ever th“ current production configuration that 
gas turbine main propulsion in uhe iifa-cycla analysis. 



a comparison is shown in Figure 1.9, also from the NIVSFA 
report [Bef« 7]. 



Additional propulsion plant arrangements are given 
in [Ref. 5] and the total number possible is only limited 
truly by the designer's imagination and cr®ativity within 
the confines of the hullform. 



3- Motor/Generatcr Los ses 



The concept of losses within the motor and the 
generarcr are comparative in magnitude and virtually iden- 
tical in nature. The losses will, therefore, be discussed as 
a group which is applicable to both the motor and the gener- 
ator. The approach of an energy balance is useful here to 
introduce these losses: 



I Energy into device] 
I on the I 

I mechanical side | 



lEnergy into devic® 
+ I on the 

I electrical side 



I Increase of | 

I stored kinetic/ | 

I potential mech. | + 

I energy in device | 



Increase of | 

stored energy | 

in the elec.&mag. | + 
fields of device | 



Heat I 
gen erated] 
within I 
device I 



( 1 . 1 ) 

The net quantity on the right-hand side of the equa- 
tion is small and positive; the majority of the energy into 
the device is transmitted through the device as an output. 
As an example, when considering rhe motor, the energy into 
the device on the mechanical side is negative in sign and 
approximately equal to the electric energy into the device 
which is positive in sign. This amount must b® balanced by 
the quantity on the righr side of the equation. 
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Figure 1.7 Integrated System Candidate, Fixed Frequency. 
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Figure 1.8 Integrated Systea Candidate, Variable Speed. 




COMPARATIVE DESTROYER SIZE AND COST* 
DD963LIKE MECHANICAL VS INTEGRATED ELECTRIC MACHINERY 
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Figure 1-9 Comparative Destroyer Size and Cost. 



Iha terms cr the right side ere all usually small ar.d may be 
zero. The iirst tera represents rhe energy associated with 
^ ch3.r.c^ 5.nd r'^^ults trorn tihs of 

Ir. the steady state it is zero. The second rerir. would result 
from a change in current to the device and the electric 
field that existed would be a source of energy (the induc- 
tion effect) . The third term is the main point of inrerest 
here. It is comprised of the following: 

* Mechanical friction losses 

* Air friction losses (windage) 

* Hysteresis and eddy-current losses 

* Resistance losses in the armature, and 

* Exciter losses (in the case of synchronous machines) . 
The object for the electrical engineer is to minimize these 
losses tc the greatest extent possible. Since it is impos- 
sible to eliminate all these losses, the problem that 
remains is the removal of the resulting heat within the 
device [Ref. 8]. The size and weight of the device is 
largely dependent on the ability to remove this heat, since 
the capacity of the device is limited by the insulation used 
in the device. The method of heat removal is not limited to 
a single mode, but is a combination of many modes of heat 
transfer including: 

* Conduction: stator to casing 

rotor to stator (through air gap) 
rctor to shaft 
shaft to casing 

shaft to external device (coupling, 
reduction gear, etc.) 

* Convection (both forced and natural) : 

casing to ambient 

stator tc internal air (or gas) 

rctor to internal air (or gas) 

shaft surface to internal air (or gas) 
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shaft surface to external air 
* Hadicticn: external surfaces to sinks 

internal surfaces zo casing 

These modes are not all very controllable. Large 
amounts cf research have been devoted to minimizing these 
thermal resistances. Efforts ro improve convecticn heat 
transrer have found than the resrrict io ns on the rcror- 
stator air gap have prevented exploitation cf this mcde as 
discussed as early as in 1926 by G.E. Luke [Ref. 9] and as 
late as in 1979 by O.N. Kostikiv, et al. [Ref. 10]. The 
improvement of the conduction heat transfer has been 
improved by minimizing the thermal resisrance cf the insu- 
lating compounds and the materials used in the construction 
cf the device, including installing thermally-conduct ive 
materials in the coil ends. 

The improvement of convection heat transfer is still 
being examined by using liquids and gases to cool the rotor 
and the stator of these devices. This procedure involves 
piping a liquid through the stator to remove heat by forced 
convecticn, and utilizing a material with high thermal 
conductivity directly adjacent to the piping to conduct it 
to the piping where it is removed by the passing liquid. 
This method is equally applicable to the rotor, provided the 
problem cf the rotating seals can be overcome in order to 
channel the liquid from an external sink to the moving rotor 
and back out. The advantage of doing this over standard 
cooling schemes is that the capacity of a device may be 
increased due to the lower internal temperatures that could 
be maintained (or the physical size and weight of the device 
could be reduced fcr an existing capacity) as shown in 
Figure 1.10 for a 40,000 hp , 180 RPM synchronous motor. The 
use cf a liquid in a bouyancy driven closed-loop, called a 
ther mos yphcn, is also possible. The amount of heat transfer 
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in this s'5thcd is dapsnder.t on 
'hs liquid being used, the x 
pi? 9 rssiso£nce in the loop, 
the heat-transfer areas then 



the specific heat capaci'y 
amper at are rise allowed, rae 

required nay be prohibitively 



large . 



Ihe two-phase cooling method within a closed loop, 
in the fora of a heat pipe, has also bean extensively inves- 
tigated. This method has an advantage over the conventional, 
lig u i d-ccol ing method in that it eliminates the necessity 
for rotating seals. These seals can be weak poinrs in the 
design from a reliability stand point, especially if high 
RPhs are involved. The heat -pipe mathod typically involves a 
confined fluid acting as a two-phase medium for heat 
transfer, exploiting the phase change as a vehicle for 
substantial heat transfer in order to remove heat from the 
rotor or the stator cr both. This method then transfers the 
heat eirher to the ambient air or to some other gas, liquid 
cr solid heat exchanger. The configuration most suitable 
for use within a motcr/gen erator on the order of 40,000 hp 
would be either to ambient air or to a fresh water heat 
exchanger external to the motor in the case of stator 
cooling, or to internal forced air through extended surfaces 
in the case of rotor cooling. 



C. PBOBIZM DEFISITICH 

The purpose of this thesis is to analyze various methods 
of cooling motors and generators of warship propulsion 
plants and to discuss the economic and reliability perspec- 
tives associated with each system. The next two chapters 
will discuss the theory and application of these two 
methods: liquid cooling of the rotor, both forced convection 
and the thermosyphon, and two-phase cooling of the rotor. 



34 



A cip p 3, !T3. S "fOT V3'^30U^ 0 ^ 1 3 rr 

coolir.g sch'rines has taen dasignsu and is currently under 
constructicn an nhe David W. Taylor Maval Ship Research and 
Developaenr Center, Annapolis, Maryland. It is sketched in 
Figure 1.11 and will be used for the analysis within this 
study. The result s/reccni aendations herein may thereby be 
experimentally evaluated. 

For the purpose of this arbitrary configuration, th«= 
model for analysis is as follows; 

1) The device will be a water-cooled frame, 
synchronous AC generator, 25,118 kVA 3 3,600 RPM. 

2) A specific conduction bar shall be analyzed for a 
typical load for all cooling configurations. It has a 
length of 0.9144n and is 0.0116m sguarewith 

a 4.763mm diameter hols bored through the length 
(herein referred to as the tube length and tube 
diameter) . The losses per bar in this configuration 
are approximately 50 w. 
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Figure 1.10 ggt£J„Co|gg|ision uith vuricus 
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II. 3ISGL^£HA52 FIG ID COOLING 

Single-phase cooling has been extensively investigated 
with regard to rotating references. Work was stinulated 
early by the printing industry, with the need to cool 
printing roll nills. The work was further encouraged by the 
gas turbine industry with the need to cool high-temperature 
turbine blades, which rotate at high RPM, and by the elec- 
trical machine industry, where the need to cool large, 
powerful generators became a necessity. Research in the 
study of flow and heat transfer in a rotating reference 
system is applicable to forced-convection flow, buoyant flow 
in thermcsy phons , and the single-phase flows of closed two- 
phase devices. 

A theoretical analysis of flow in a heated, horizontal 
tube without rotation by Morton [Ref. 11] was completed in 
1957. Although his analysis was limited to non-rotating 
flows, it gave considerable insight into buoyancy-induced 
secondary flows. His analysis was limited to laminar 

convection in horizontal, uniformly-heated tubes at low 
Rayleigh numbers (based on the temperature gradient along 
the pipe wall) . The work was an extension of the study by 
Nusselt who ignored the gravitational (acceleraticnal) 
effects and, therefore, was independent of orientation. 
This study was a numerical solution to the governing equa- 
tions assuming small rates of heating and a notable varia- 
tion of density along the pipe length giving rise to 
secondary flows. A similar study was completed by Kays 
[Ref. 12] in 1955. Morton studied this secondary flow 

driven by the buoyancy effects and the resultant temperature 
distribution by successive approximation in a power series 
(this is similar to the flow in a pipe rotating about a 
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perpendicular axis). Secend-ordar terms ware required to 
account for the secendary flow. The basic aquations and 
davo Icp rent art prtsant.sd in Appendix k, and are the basis 
for many later studies. This development is adapt<=d to th<= 
rotating case by the substitution of centripetal accelera- 
tion for the gravitational acceleration term. 

The configuration of the cooling loop has three distinct 
regions which must be considered; the radi al sections of 
tubing (oriented along the line of the centripetal accelera- 
tion) , the ent ran ce regions of the axial sections, and the 
axial sections through the conductor bars (oriented parallel 
to the axis of rotation). Additionally, it must be deter- 
mined whether the flow is 1 a aina r or t urbulen t. The devel- 
opment will be reviewed chronologically within each topic 
with sample calculations made for pertinent correlations in 
Appendices E and C as noted. 

A. HISTORICAL DEVELOPMENT OF ONCE-THROOGH COOLING 

In the 1968 publication, "Recent Advances in Heat 
Transfer" [Ref. 13], Kreith made a survey of virtually all 
the work to that date that had transpired in the field of 
convection heat transfer in rotating systems. His section 
cn concentric cylinders and rotating tubes summarized the 
works mentioned herein, although principally evaluating the 
reports using gases as the fluid media. He presents the 
concept of utilizing a thermosyphon for the cooling of 
rotors, the theory of which is intimately related to the 
previous work with single-phase heat transfer, which will be 
discussed later. Kreith also describes the three different 
sections (radial, entrance, axial) of the device which this 
study is addressing. These have been analyzed individually 
both expermentally and theoretically. The inlet region has 
the swirl effect due to the Coriolis force, which is a 
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function of t.he Rosscy nuaber, or th? swirl number (rh = 
reciprocal of rha Eossby number) ; tha radial and axial tube 
sscticns have, whan a significant centrifugal accel<=ratio'' 
term is present, noticeable hydrostatic effects and, addi- 
tionally, buoyant effects whan a density differential 
exists. 



Extending his research to a vertical tubs rotating 
about a perpendicular axis, Morris [Ref. 14] confirmed that 
the approximate solutions must include second-order terms by 
studying low heating rates and using a method similar to 
Korton, but including tangential terms as well as radial 
acceleration terms. He identified the three components of 
convection: 

1) Forced convection due to the externally-generated 
pressure gradient, 

2) Gravitational buoyancy in the axial direction, and 

3) Rotational buoyancy due to centrifugal and Coriolis 
forces. 

Ihe relative magnitude of the terms due to item 3) is char- 
acterized by the Rossby number. This is a ratio of th<= 
inertia force to the Coriolis force. Although the Rossby 
number exists in the velocity and temperature fields for 
solutions through the second order, it has no effect on the 
solutions for heat transfer or flow resistance. Since at 
the nominal rotation rates and eccentricity that exist in 
rotating electric machinery, the centrifugal acceleration is 
much much greater than the gravitational acceleration, the 
former will dominate. 

In their 1968 report on heat transfer in rotating, 
radial, circular pipes in the laminar region, Mori and 
Hahayama [Ref. 15 ] continued their research by evaluating 
steady, fully-developed flow with analytical techniques. 
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Usir.G *he toun dary- layer ccr.capt, they correlated ^h=ir 
results vi+h ex peri niental work by Trefether. They gav* 
corr elauicr.s for ur.s (du/liuo) rario as a funcuior. of 'hr 
argular velocity, Reynolds, Prandtl, and rotational Reynolds 
numbers for fluids with Prandtl numbers above and below 
unity (i.e., gases, liquids, glycerol, and liquid metals). 
Their snudy confirms the stabilizing effect of the secondary 
flow cn the critical Reynolds number and notes that laminar 
solutions are still valid when the Reynolds number is 
considerably large (12,600 3 Q. =500 rad/s). 

As a second report on the subject of rotating radial 
pipes, Mori, et al. [Bef. 16] evaluated secondary flow 
effects due to the Coriolis force by assuming a boundary 
layer along the wall. with experimental results and corre- 
lations from previous work, they were able to show that the 
effect of secondary flow is less in the turbulent region 
than in the laminar region, although still on the order of a 
lO’i increase over the value without secondary flow. The 
correlations they presented agree with their previous report 
[Bef. 15] and with experimental results with air; they are 
presented as a ratio of Nusselt numbers based on an exponen- 
tial function of the Reynolds and Prandtl numbers and are 
given for a wide range of Prandtl numbers including gases 
and liquids. 

Ito and Nanbu [Bef. 17] experimentally measured 
friction factors for flow in rotating, straight pipes of 
circular cross-sect icn, and compared their results to those 
of Trefethen and Mori 5 Makayama. Their empirical correla- 
tions for friction factor ratios for two ranges of Kt (where 
Kt = J2/5e and J= d2/ v ) are shown in equations (2.1) 
and (2.2) for the turbulent region, Kt > Ktcr. These are as 
follows : 
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0.942 Kto.282, 1< Kt <500, 



( 2 . 1 ) 



\ 0.942 Kto-05, Kt>500. (2.2) 

Eelow a Kt of 1, the fricticn factor practically ccincidas 
with th .2 stationary friction factor givan by Blasius: 

f ^ = 0.316 Re-0.25 . (2.3) 

For laminar flow, which occurs when Kt is below Ktcr, the 
following correlation applies: 

f/f = 0.0883 Klo*2S(1 + 11.2 Kl 
^1 

where Kl = J*F.e. Equation 
2.2*102< Kl <107 and (J/Re)<0.5. 
without rotation is given by: 

f . = 64/F,e (2.5) 

^1 

Transition occurs at: 

Ktcr = 1. 07 J1.23* 1C3 (2.6) 



-0.325) , 



(2.4) 



(2.4) is valid for 
The friction factor 



or 



Re =1.07 JO . 2 3«i 03 (2.7) 

cr 

for the range 28< J <2*103. Calculation of the pressure 
drop would be required in the detailed calculations for 
design of a cooling system, especially of the closed, 
rot at ing-loop system. 



2 . 



Axial SactioTiS 



Kuo, 3 t al. [Haf. 18] analysed ths hoa t-trar.sfar 
characteristics of water flowing through a partially- filled 
pipe which was rotating about its own axis and which was 
maintained at a constant heat flux. They also summarixsd 
the previous pertinent works that wera available at the time 
the paper was written. The paper illustrates hew the tran- 
sition from gravitational to cenxrifugal force is accompa- 
nied by instabilities as previously discussed by Taylor 
[Ref. 19] and [Ref. 20]. The authors calculated a critical 
Taylcr number at which the Nusselt number changes its func- 
tional dependence on the Reynolds number. The paper also 
shows how the heat transfer is accomplished within a layer 
around the periphery of the pipe. A correlation for this 
film thickness is provided, and experimental results confirm 
its validity and illustrate how the increasing depth of the 
film in the rotating pipe, up to completely full, does not 
enhance the heat transfer beyond the critical depth value. 
The range of the rotational Reynolds numbers used in the 
experimental work was between 4,000 and 20,000; the enhance- 
ment (Nu/Nu^) was approximately up to 240?. 

Oliver [Ref. 21] analyzed natural convection in a 
horizontal tube with no rotation, following the work of 
Colburn [Ref. 22]. This study, while pointing out several 
problems with the original research and basic criticisms 
noted by Martinelli [Ref. 23], combines the work of many 
authors to obtain a correlation (equation (2.8)) for the 
Nusselt number in terms of the Graetz, Grashof, and Prandtl 
numbers and the L/D ratio. The paper also points out that 
the Grashof number must be evaluated using the temperature 
gradient along the pipe axis instead of the gradient from 
the pipe wall to the fluid mean temperature. 

NU 1.75 [Gz^+5.6*10--»* 



m 
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( 2 , S( 



(Gr ?r (L/d) 5 0 . 0 . 3 33 . 

m m 

"Ihe correlation is supposed to incorporate the sfiecrs of 
natural convection along with forced convection by vector 
addition, but he nctes that the velocity profile of the 
fluid flew is critical in the analysis and that natural 
convection in the case of horizontal tubes should be inde- 
pendent of L/D beyond the entry region. For this concept, a 
correlation is also provided for th'=' Nusselt number withou- 
the L/D term: 

Nu (yyyj^) 0 .i 4 = ■'• 75 {Gz^ 

+ 0.C083(Gr Pr )o. 7550 . 333. (2.9) 

m m 

The study of the Taylor vortices was continued by 
Pattenden [Ref. 24] using a horizontal (axial) tube rotating 
about its own axis. Due to problems in the slip-ring appa- 
ratus used in the experimental work, the errors in the accu- 
rate measurement of the temperature precluded specific 
determination of exponents for the Reynolds number in the 
correlation for the hea t-transfer -coefficient equation. 

Qualitative analysis of the enhancement due to tube rotation 
was made, noting that the axial-flow effects were small 
compared to the rotational effects. Taylor vortices were 
neither confirmed, nor denied; the flow was sufficiently 
turbulent to mask the effects that the vortices may have 
had. Rotation rat®s were varied to a maximum of 4,000 RPM. 

Mori, et al. [Ref. 25] investigated further into 
the effects of buoyancy on forced-convective heat transfer 
in horizontally-oriented, heated tubes with laminar flow. 
His study extended the work of Graetz and Nusselt, and went 
beyond the limitation of small heat fluxes such as those 
srudied by Colburn, McAdams, Martinelli, et al., Jackson, et 
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ai . , and Oliver. The sajcrity of the experimental vorl< bv 

the authors was performed with a constanr wall tetrperat ’ire 

and large L/D ratio. Their study involved a large neapsra“ 

ture difference between the wall and the fluid such as would 

occur with large heat flux. The change in viscosity with 

tenperanure is appro ximanely covered with the (y /y )o.i4 
. . * " w b 

tern rn equation (2.8), but the effect of natural convection 

has not been included specifically. This effect was shown 
[Ref, 11 ] as a function of the ReRa product; it is appli- 
cable in the range of ReRa less than 3,000, but within this 

range the effect of natural convection does not exceed 

several percent. The Rayleigh number is a function of the 

fourth power of the diameter; therefore, the effect that 
natural convection makes in small-diameter channels is 
small. This illustrates that in the thin layer near the 
wall, the temperature and velocity change abruptly and are 
quite different from that in the inner part of the fluid. 
When the ReRa product is large, the secondary flow is strong 
and the velocity distribution in the axial direction is 
urterly different frcm the Poiseuille flow assumed in the 
previous works. The Nusselt number is then calculated by; 

Nu = 2aq/ k (t -t ) (2.10) 

w m 

where the stationary Nusselt number for constant heat-^’flux 
ca s e is; 



Nu g = 48/11. (2. 11) 

The experiments were conducted with air and fluids with 
higher Prandtl numbers; they were only considering heating 
in a gravitational field (not an acceleration field due to 
rotation). In this case, the layer near the wall becomes 
thinner and the Nusselt number increases with the increasing 
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HeRa. Ths local Nusselt nuitib<;r for the experimental ccrre- 
laticn (for air) is given by: 

Nu = 0.6 1 {?.8Ha)o. 2 (1+1.8/ (RaEa) 0.2) (2,12) 

This study vas extended into the turbulent region, where the 
highest point of the temperature distribution was found to 
shift in the directicn of gravity slightly. In this case, 
the effect cf the secondary flow was not as significant as 
in the case of laminar flow. The Nusselt numbers calculated 
within this regime were in good agreement with those of 
Colburn [Ref. 22], where: 



Nu = 0.020U Reo.a. 



(2. 13) 



Mori and Nakayama used the theoretical analysis of 
their previous research [Ref. 26] and [Ref. 27] on straight 
pipes rotating about a parallel axis [Ref. 28], Analyzing 
the body forces driving the secondary flow caused by density 
differences in the centrifugal field and the Coriolis force 
with regard to the flew, they used fundamental principles to 
characterize the flow and temperature fields. This was done 
assuming an effective secondary flow due to buoyancy and 
using the rotational Reynolds number. This would indicate a 
rapid divergence frem the results of Morris [Ref. 14] and 
shows that Coriolis effects cannot be ignored. The paper 
gave correction factors as a function of Reynolds, rota- 
tional Reynolds, and rotational Rayleigh numbers, and 
concluded that the Coriolis effects diminish with increasing 
eccentricity (higher g- fields) to a negligible value when 
the ReRa product is greater than 10,000,000. 

In 1968, Nakayama [ Bef . 29] further analyzed a hori- 
zontal, straight pipe rotating about a parallel axis. He 
assumed, as previously [Ref. 15], that an effective 
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secondary flow in fully-developed conditions of flow and 
teniperature fields existed, and included body-force terns 
explicitly. The correlations he obtained for both the rario 
cf friction factors and the ratio of Nusselt numbers are 
valid over a large range of Prandtl numbers for both liquids 
and gases, as well as for a wide range of Reynolds and 

Grashcf numbers. The correlation for the Nusselt number for 
liquids in fully-developed flow is: 

Nu = Rso-9 Pro-'* {0.033 ( Re/ r ^ • s) 1/3 o 

(1 + 0.0m/(Re/ r^-®) (2.14) 

where F = Re 22/1 3 (GrrPr® • 2/1 3 . p term is the 

ratio of inernia force to body force. Inertia force is 
represented by Re22/i3^ and the remaining terms in F 
represent the body force. The body force is either the 

Coriclis force in the case of the inlet region, curved 

regions, or the radial arms, or the centrifugal force in the 
case cf the axial section. The numerical result of his 
correlations for the model presented in Chapter I is 
included in Appendix B for both turbulent and laminar 
liquids. 

Sieawarth, et al. {Ref. 30] further evaluated the 
effect rhat bou ya nc y-induc ed secondary flow would have on 
laminar flow in a heated, horizontal pipe without rotation. 
His work shewed clearly that secondary flows did exist due 

to density variations and that the heat transfer was 

enhanced by them as a function of the Grashof and Prandtl 
numbers. This affect is somewhat negated in the presence of 
extreme acceleration as that in the periphery of a rotor, 
but remains nonetheless, provided the turbulence is below 
the critical Reynolds value. 
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!*ith no assumptions as to the structure o£ the flew 
and temperature fields to simplify the governing equations. 



W Sm/ -i W <-4 



Ilcrris oresented a 



of the rotating cylindrical tube [Ref. 31] with laminar, 
fully-developed flow. Data were compared to the theoretical 
results for the case of air, water, and glycerol. The 
correlation obtained is in terms of the 

Eayleigh-Reynold s-Pr andt 1 product as is given by: 



Nu = NUg {0.262 (EaRePr) 0- 1 ^3} . 



(2.15) 



The value for the test model is calculated and presented in 
Appendix B. 

Stephenson again studied this parallel, rotating- 
pipe, heat-transfer problem for fully-developed turbulent 
flow. He compared his experimental results and correlation 
with the earlier work of horris and Woods [Ref. 32] and they 
compared favorably within the turbulent entrance region. He 
noted that the rotational buoyancy was not a strong factor 
in the secondary- flow in the fully-developed region. 

Because of this, his results, when compared against 
Nakayama's results assuming a strong effect [Ref. 29], did 
not compare well. His correlation is listed below: 



NU = 0.0071 Re0*®8 JO. 023 



(2. 16) 



and is compared with previous correlations for the test 
model in Appendix B. 

3 . Ent rance Regi ons 

In 1969, the Institution of Mechanical Engineers, 
Thermodynamics and Fluid Mechanics Group, sponsored a sympo- 
sium, on the subject of heat transfer and fluid flow in 
electrical machines. Included in the presentations was a 
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t ate r-rotor test rig for avaluating tha hsa t- tr-insf 9 r 



that enhancement in the turbulent flow regimes was net as 



■•■he reason was the inability of the hydraulic-diameter 
concept to accurately account for the flow ccndifior.s in the 



turbulent correlation with theoretical prediction when the 
cooling scheme was modified by shortening the axial path. 
This effect, although not explained, could have been due to 
the Coriolis swirl effect in the inlet region enhancing the 
heat transfer above that to be obtained in the fully- 
develcped region. 



Prandtl numbers such as those of gases, are very difficult 
to eliminate and the Corio lis- effects on the heat transfer 



Eeynclds-Rayleigh product increases above about 1,000,000, 
there is a tendency for the amount of enhancement to 
diminish. This is attributed to the turbulent effects over- 



ticn on heat transfer is to enhance it in the entrance 
regions by the secondary flow due to the Coriolis effects 
and in the fully-developed, horizontal-tube regions by buoy- 



by Petukhov and Polyakov in the USSR with much the same 



entrance region of a tube, where the Coriolis terms domi- 
nate, is addressed by Morris and Woods [Ref. 32]. 
Correlations are presented for air and are applicable for 



laminar and the turbulent cases, are presented as functions 



characxeristics of hydrogen cooling 



His testing shewed 



high as theoretically predicted. It was hypothesized that 



non-circular ducts. The tests were able to give a good 



The effects of the entry length, especially with low 



may be noted even at high Reynolds numbers 



As th e 



iding the secondary flow enhancement. The effect of rota- 



ancy effects. A similar survey of the technology was made 



results [Ref. 34]. The heat-transfer problem in the 



ether gases as well. These correlations, for both the 



of the product 



Reynol ds 



and the 



4 9 



of th e 



num ber 



rctational-F.eynolds liumber . Th9 Morris and '■icods oarar 
notes that further work is required in the case of liquids, 
although a siriiiar approach is valid. Mori and Uakayanic. 
also presented a survey of the stat a- of- the-ar t technology 
regarding tha heat-transfer characteristics of rotating 
pipes and ducts [Ref. 35], The questions of the entrance 
effect and its length were emphasized as bair.g generally 
unanswered, but they gave a general guideline for deter- 
mining tha extent of the entrance region as approximately 20 
times the diameter of the pipe. They also verified correla- 
tions of previous works for helium and water. 

4. Com bin a t io ns of Ra di al. H oriz ontal , and Ent rance 

Sec tio ns 

Also presented in the 1969 Institution of Mechanical 
Engineers Symposium was a study by Lambrecht [Ref. 36], 
which discussed the problems associated with water cooling 
of rotors. His paper summarized the theoretical considera- 
tions, as well as reviewed his previous work and the work of 
fellow German researchers (Neidhoeffer and Ingenieure) . He 
noted the superior cooling properties of water and the size 
and weight reduction, along with the improved efficiency of 
water-cooled machines. He listed values for the optimum 
cooling duct size based cn the heat-transfer and 
electrical-loss characteristics for various gases. The 
method used for optimum duct calculation was based on the 
principle that the pressure drop in the duct decreases with 
increasing duct size and electrical losses increase with 
decreasing conductor area (increasing duct size) . A suit- 
able calculation of this type is necessary in any final 
design fcr cooling of these devices by any method devised 
for either the single- or the two-phase approach. 

In 1970, Sakamoto and Fukui measured heat- and 
mass- transfer coefficients for air and oil, specifically for 
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use in ceding rotating electric machines (Bef. 37]. Noting 
that the most effective and economic method to increase the 
*3 C V ^ IT C P C — C Z* I- C 3, 2. rn 3. C h. ^ n 3 S 3 S ^ O ^ ID C ^ JTt 

cooling of the insularing marerial (especially in th=> rotor 
conductors and drums) , a geometry similar to a typical rotor 
was fabricated and suitable coolants were tested. The 
convective heat transfer was characterized as a function of 
the L/d, Re, 3rr, Pr variables, with the H/d parameter held 
constant. When the shape factor of eccentricity H/d is 
sufficiently large, the centrifugal acceleration may be 
assumed to be a uniform acceleration { 2H) acting on the 

tuba axis, which is similar to the force gravity exerts on a 
horizontal tube. with a lew Reynolds number, the empirical 
correlation was in the form: 

\ 

Nu = NUg (1+C (Ear) ^/Gz) (2.17) 

where the values of the constants were: 

C = 0.03 m = 0.75 n = 1/3 

and where Nu = 1 . 8 6 (Re Prd/L) o • ^ 3 for large L/d values, 
and Nu = 2.67 gzo* 33 for small L/d values. 

This was valid within +10^ to -17% for Reynolds numbers 
between 162 and 2,700 and rotational Rayleigh numbers to 
20,000,000. It was also found that the length of the 
entrance region seemed to have no effect on the correlation, 
which wculd imply that a secondary flow existed in both 
configurations. An interesting point in this paper is that 
the rotational Grashef number (hence rhe rotational Rayleigh 
number) was calculated using the temperature difference 
between the wall and the fluid instead of the axial tempera- 
ture gradient. 

Woods and Morris investigated laminar flow in the 
rotor windings of directly -cooled, electrical machines in 
1974 [Ref. 38] and noted that if water was used as the 
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coolaTit in li^u of hy3roc'?n (which is currently widsly 
used) , an efficiency improvement of >0.5/5 could be o'ctainsd. 
Given that the 0.5"( cf a large, say 1,000 KW, generator was 
the only consideration made (and neglecring the size and 
weight reduction and the extended operating life possible 
due to the cooler operation) , they proceeded to analyze the 
fundamentals of the problem. They surveyed the previous 
works directly applicable to this area of coding; 

1) Morris [Ref. 14] used a series-expansion technique 
(valid fcr low rotation rate and low heating rates) , and 

2) Mori and Nakayama [ Hef . 28] assumed a secondary flow 
(which was claimed to be valid fcr high rotation rates) , and 
used an integral-type analysis. 

Hoods and Morris claimed both were inadequate due to the 
restrictions imposed by the nature of their solutions. They 
attempted a versatile and "exact" solution for laminar flow 
in the fully-developed region by solving rhe governing equa- 
tions wi-^h a numerical procedure. Numerical solutions were 
presented for both the friction factors and the Nusselt 
numbers. Their analytical values were compared with experi- 
mental results and to correlations of previous works. 
Discrepancies were explained and the difficulty in obtaining 
fully-developed conditions in the experiments of their work 
and previous works was emphasized. The axial density- 

variation affects were noted, also, as a potential for error 
in their analysis at high values of rotational Rayleigh 
numbers. 

Nakayama and Fuz icka evaluated the generator 
problem, where reasonably large radii of rotors (approxi- 
mately 1 m) and high RPM (3600) create high centrifugal 
acceleration affects [Ref. 39]. This acceleration causes 
the centrifugal buoyancy term to be significant. It also 
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esph'3 2iz=s the Coriolis term i 
curved inlet regions. They 
ratio of friction factors and 



r. tha radial tubs segssn t~ and 
presented correlations for the 
the ratio of Nusselt runbers. 



arriving at friction-factor rarios for each of the three 
regions (radial, curved inlet, and axial) , and compared the 
results to experimental data. Their correlations are listed 



below : 



f/fg =2.2 Ro“ 0*33 (radial) 



(2. 18) 



f/fo = 1.5 Ro“0*3 (coil end) 



(2. 19) 



f/fo = 14.0 Ro—3. e (axial) 



( 2 . 20 ) 



Nu = ReO-3 Pro-* {0.014(Re/Ro2-5) 0.124) (radial) (2.21) 

Their report confirms the results of Nakayama's correlation 
(equation (2.14)) for the axial sections. 

The application of the technology of rotating 
cooling schemes is net limited to the printing industry and 
electric machine area; extensive research has been devoted 
to this field of study by the gas turbine industry. The 
Advisory Group for Aerospace Research and Development 
(AGARC) held a meeting in September, 1977, at Ankara, 
Turkey, to discuss High Temperature Problems in Gas Turbine 
Engines. Numerous papers were presented on the subject 

including a paper by W.D. Morris on flow and heat transfer 
in rotating coolant channels [Bef. 40]- He used a selection 
of experimental results to illustrate the influence of 
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rotation cn h<=at transf^ar, and daraonstrated that Coriolis 
and csntripetal inartiai effects significantly altar ths 
heat -transfer characteristics ccnparad to the stationary 
case. Using the governing equations, he illustrated how the 
Navi er-Stokes equations are modified to include the acceler- 
ation effects of the Coriolis forces and the centripetal 
forces. For constant properties, these equations then show 
hew at distances sufficiently downstream from the =ntry 
region, where the axial gradients of the velocity are negli- 
gible, the eliminaticn of the pressure-gradient terms from 
the radial and tangential momentum equations causes the 
Coriolis terms to vanish identically as a source for the 
creation of secondary cross-stream flow; see also [Ref. 14]. 

In the entry region, the Coriolis terms interact 
with the developing axial velocity and create secondary 
flows perpendicular tc the axial direction, even with the 
constant property condition. Finally, in the fully- 

developed region, the effect of the rotation manifests 
itself only as a cross-stream pressure distribution if buoy- 
ancy terms are neglected. If the density is allowed to vary 
with temperature, the centripetal acceleration terms need 
cnly be included as the Coriolis terms do not affect this 
action. This cross- stream, buoyancy-induced flow gives rise 
tc greatly-enhanced beat transfer and an attendant increase 
in flow resistance. The effects of the eccentricity, with 
regard tc the Coriolis terms, has little effect, provided 
that the ratio of the radius (H) to the diameter of the tube 
(d) is greater than 5 (only very small radius of rotation 
will have any noticeable effects). 

Marto [Ref. 41] lists the papers presented at the 
14th Symposium of the International Centre for Heat and Mass 
Transfer (ICHMT) held in Dubrovnik, Yugoslavia, in 1982. 
Papers of particular interest to this study included a paper 
by Jehnsen and Morris [Hef. 42], wherein the concept of 
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csntripetal-acceleration-induced secondary flow was 
to be cf little importance, and its influence was s 






be totally hydrostatic. 7hair work also confirmed the 
effect cf Coriolis acceleration with regard tc secondary 
flows, especially in sections with relatively small L/d 
ratios. No quantitative recommendations were presented. 



5 • rransit ions 



Mori [Ref. 25] states, when considering turbulent 
flow, that it is not necessary to consider the influence of 
buoyancy on heat transfer. The secondary flow was found to 
suppress the turbulence level when the turbulence at the 
inlet region was high and an empirical formula for rhe crit- 
ical Reynolds number, in terms cf the Reynolds number and 
Rayleigh number product, was calculated: 



Re = 123(ReRa) o- 2s. 
cr 



( 2 . 22 ) 



Contrary to this, when the turbulence at the inlet region 
was low, the critical Reynolds number was higher (7,700 vs. 
2,000), and heating decreased the critcal value. The net 
result was that when the ReRa value was high, the secondary 
flow caused by buoyancy makes the critical value of Reynolds 
number tend toward the same value, whether the turbulence in 
the inlet is high or low. In the case of low turbulence in 
the inlet, the value for the critical Reynolds number was 
given by: 



Re = Re /(1+0.14 ReRa*10-5) . (2,23) 
cr cr, 

The graph shown as Figure 6 2 in the Krieth article [Ref, 13] 
is shewn below as Figure 2,1. This graph shows the heat- 
transfer as a function of both the RPM and the flow-through 
Reynolds number. It is easy to see the flow's transitional 
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dep'snd^nce 
ths oxpsri 
Har=7,40G, 



on th9 Reynolds number as a function 
mental model being evaluated, where 
the transition is delayed to Re=11 



of PPM. 
Re =7, UOO 
000 for 



cor r elan ion 



that assumes high turbulence at the entrance 
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Figure 2.1 Effect <?f Rotational Speed 
in a Rotating Pipe. 



The transiticn to turbulent flow is suppressed by 
rotation and has been g ualit ati vely noted by numerous 
authors. An exact correlation for transition to turbulent 
flow in a r ctating“reference frame with heating is net known 
to exist. From the available literature, it is clear that 
the transiticn follows a path from the laminar correlations 
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to the turbulent correlations. This path, in the cas^- of 
the rotating reference, has a ’'dip" as shcvn qualitat ivelr 
in a plot of the Nusselt number against the Reynolds numb^^-r. 
Figure 2.2, whereas in the non-rotating case, no "dip" 
occurs. 




Rej_ > Re2 > Re^ > RS4 

Figure 2.2 Transitional Characteristics of 
Hotatxng Systeus. 

6. Summary 

Eased upon the information reviewed in the litera- 
ture, it may be concluded that the three regions within the 
motor/gener ator application (shown in Figure 1.11) have very 
different heat-transf er characteristics. 

1) In the radial section cf the coolant path, the Coriolis 
forces tend to dominate the heat transfer and friction. 
This effect is to enhance the heat transfer and increase the 
pressure drop in these sections. The heat-transfer enhance- 
ment is substantial and could be effectively exploined for 
the heat rejection from the device by axial forced-air 
convection of the substanti ally-hollow, synchronous rotor. 
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2) Ir. ths coii-sr.d regions (bcth at tht radial conr. ioa 
points and tha coil-end loops) the Coriolis force is strcna 
and is coupled with the strong, cantri petal forces. rhs 
result is enhanced heat-transfer and friction-factor 
effects, but to a lesser degree than in the radial sections. 

3) In the winding bars, the centripetal forces doicinate: 
these forces result in hydrostatic effects which limit flows 
that would be affected by the Coriolis forces and result in 
some degree of seconcary flew due to buoyancy. The amount 
of secondary flow is questionable and is somewhat affected 
by the degree of turbulence and the amount enhancement is 
less specific. 

The enhancement in the laminar-flow heat transfer 
has been reported by numerous authors. This can be seen in 
Figure 2.3. 3oth the Nakayama [Ref. 29] and the 
Scods-Morris [Ref. 31] correlations for Nusselt number fall 
between the classical laminar value of U8/11 and the turbu- 
lent value of the D ittus-Bo elter correlation. They have the 
same trend rhroughout the range of RPM shown and differ by 
only a few percentage points. Apparently, the effect of 
rotation in laminar flow enhances the heat transfer because 
of intense secondary flows due to centripetal forces. An 

interesting point is that, at the higher RPM values, the 
correlations of Nakayama, Woods-Morris , Stephenson, and 
Dittus-Bcet ler all are exceptionally close together. 

The heat-transfer predictions for turbulent, fully- 
developed, parallel flows are shown in Figure 2.4. This 
figure illustrates that the resultant heat transfer is some- 
where in the vicinity of the classical turbulent correlation 
of Di ttus-Bcel te r for the current model. 
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Figure 2.3 Variation Qf. Latainar Nusselt 
Correlations with RPH. 
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RPM 

Fiqure 2.4 Variation of Turbulent Nusselt 
^ Correlations with RPM. 



The hea^. transfer is either scmewhat enhanced (on the order 
of 30^?) in the lower HPJi ranges (Nakayama ( Eef . 29 ] } , or 
somewhat suppressed hy the high centripetal accc It ration 
forces as shown by the correlation of Stephenson (Sef. 43 j. 

The enhancement of the heat transfer in the Nakayama 
correlation [Ref. 29] in the lower RPM range is assumed to 
be due to the residual Coriolis effects from the entrance 
regions remaining longer in the lower centripstal accelera- 
tion field. 

E. HISTORICAL DE7EL0FHENT OF THERHOSTPHON COOLING 

The thermosyphon is a device that creates the motion of 
the working fluid solely by buoyancy forces. It has a 
region where heat is supplied to the working fluid and a 
region where heat is rejected. The less buoyant, cooler 
portion cf the fluid moves in the direction of the accelera- 
tion vector, while the opposite is occurring for the warmer 
portion. The ther mcsyphon also has a third region, the 
coupling region, with problems of fluid shear, mixing and 
entrainment. Figure 2.5 illustrates the conceptual arrange- 
ments of the three basic thermosyphon systems. 

Kreith [Ref, 13 ] mentioned in his review that the ther- 
losyphon (i.e., a liquid -filled container) had a great 
potential for cooling of electric rotating machinery. In 
motors and generators, the use of a thermosyphon would 
preclude the necessity of using rotating seals, which have 
marginal reliability in applications of high RPM and volu- 
metric flow. In 1S73, Japikse [Ref, 44] gave a thorough 
review of the literature regarding thermosyphon technology. 
Beginning with a reference to the Davies and Morris paper 
[Ref. 45], which will be discussed in detail later, Japikse 
categorizes thermosyphons as follows: 
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1) the 



nature of bound a 



(either 



opener or 



mass flow) , 



2) rhe regime of heat transfer (purely natural ccr.var- 
tion or mixed natural and forced convection). 



3) the number and types of phases present, and 

h) the nature of the body force present (gravitational 
or rotational) . 

These broad categories are then restriced to "....a 
prescribed circulating fluid system driven by thermal buoy- 
ancy forces.'* 

The use of the closed thermosyphon in rotating refer- 
ences, such as in gas turbine cooling, has been studied in 
detail. In 1965, Eayley and Lock [Ref. 46] experimentally 
investigated the performance of closed thermosyphons. They 
utilired the theoretical development of the open thermosy- 
phon, modified the theory, and experimentally verified their 
results. It vas shown that the critical operating parame- 
ters include: 



1) Length-to-Diameter Ratio. This ratio controls the 
characteristics of the flow and determines the extent of the 
coupling region. 

2) Heated-Length-to-Cooled-Length Ratio. As this ratio 

tends toward zero, the analysis is identical to the open 
ther mcsyphon correlations. This restricts the use of the 

analysis for the closed thermosyphon to Lh/Lc ratio to be 
not much greater than one. 

3) Coupling Region. This is the region that compli- 
cates the analysis and causes the Prandtl number to affect 
the heat transfer. It exists in three distinct modes; 
conduction, convecticn, and mixing. These are regarded as 
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id5s.ls 3.nd actual coupling usually consists of a combination 
of all thre=’. 

Therniosyphon technology was also applied to trn nsf comer 
cooling, nuclear-reactor cooling, heat-exchanger fins, home 
and industrial furnaces, cryogenic cool-down apparatus, 
steam tubes for bakers’ ovens, internal combustion engine 
cooling, and env ircnmenta 1 control in the space vehicle 
program, as well as many others mentioned in the Japikse 
arti cle. 

In 1971, Bayley and Martin [Ref, 47] also reviewed the 
stat e-cf-the-a r t technology, with particular emphasis on 
gas-turbine applications. They studied both the open and 
closed thermosyphon systems. The use of an open thermosy- 
phon system in a rotating reference, such as in gas turbine 
cooling, as well as in the cooling of electric, rotating- 
machinpry area, has problems regarding the fluid selection 
and its saturation temperature and pressure. In the open, 
rotating reference system, the region that receives the heat 
is under high acceleration forces with high resultant pres- 
sures. As the natural flow due to buoyancy occurs, the 
heated fluid moves frcm a high pressure to a relatively low 
pressure area where a phase change to vapor may occur, 
possibly even explosively, blocking the flow. They intro- 
duced the idea of using two-phase cooling within the closed 
tharmosyphon as a way of exploiting this phenomenon, by 
using the high heat fluxes associated with the latent heats 
of evaporation and condensaticn and the lower temperature 
gradients associated with the phase change processes. This 
also has an added benefit of reduced weight over the 
single-phase liquid system. They noted that further 
research in the area of the Coriolis force effects on the 
heat-transfer chact eristics in a rotating reference needed 
to be accomplished. 



63 



Acceleration II, Field 
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The Open Thermosyphon The Closed Thermosyphon The Closed Loop Thermosyphon 

Figure 2.5 Various Thermosyphon Arrangements. 



Applyinc the technology of the r her cnosyphon tc alecrric 
oachines, Hcrris and Davies [Hef. 45] adapted the principles 
of the closed t her lacsyphon to a closed rotating loop that 
could be placed into the rotor assembly of an electrical 
device. This configuration avoids the entry choking and 
nixing that occurs in the typical closed ther :ncsyphcn and 
generally avoids the complex mid-tube exchange process, as 
well as the adverse core-boundary layer interaction. The 
strong acceleration force present in the periphery of the 
rotor of electric rotating machinery, coupled with the 
heating due to the electrical losses, induces stong loop 
circulation. This heat can then be rejected either to the 
shaft of the rotor acting as a heat exchanger itself, tc the 
air passing through the hollow rotor region and the air 
passing through the end-bell region of the rotor ends, or to 
both. The shaft, acting as a heat-exchanger itself, can be 
cooled with an internal, rotating heat-pipe (or two-phase 
thermosyphon) which could transfer the heat axially to 
another heat sink outside of the motor/generator casing. As 
an alternative, the shaft cculd be internally coded by 
forced convection of water (or some other fluid) to an 
external heat exchanger through the non-coupling end of the 
shaft by the use of an axial rotating seal. The reliability 
and the casualty -control aspects of the axial rotating seal 
make it preferable tc the radial rotating seal mentioned 
earlier. 

Theoretical analysis of the loop-flow has used simple, 
or.e-dimensicnal force and energy balances, and need net be 
repeated here. The closed loop is shown in Figure 2.6 as 

conceptually applied to a rotor. Of course, numerous loops 
would be designed into the rotor to carry away the heat 
generated. Incorporation of a one-way valve was suggested 
as being necessary tc insure one-directional flow created by 
the mean temperature difference and the pressure differences 
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Figure 2.6 Rotating Loop Theraosyphon for a Rotor. 



in the limts of the circuit. Experimental testing was 
conducted in a test rig similar to the rig proposed for the 
evaluation of the systems studied herein and now under 
construction at DTNSBDC. A modest range of rotation rates 
was employed in the Morris and Davies test rig, typical of 
most motor applications. The details of their test section 
were very similar to the model used for analysis herein; 
25. U mm (one inch) diameter copper rod through which a 6.35 
mm (0.25 inch) hole was bored. The experimental results 
were placed in the form; 
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Nu = f <Grr,Ec, Pr , Ac) . 



{2. 2U) 



The ratio of angular acceleration to gra virationa 1 accelera- 
tion (Ac) was the authors’ way of including the rotational 
effects (previous authors have used the rotational Grashof 
number) . The 
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cf cooling, if sel€ctsd, would be inappropriate for the 
actor device (i.e,, lew ?.?!!); on the other hand, it could 
he very reasonable for rhe generator. The friction factor 
was assumed to be 64/Re and no correction was made for rhe 
rotation in the friction- factor calculation; the minor 
losses were neglected. If the friction factor is increased 
due tc the rotation, as in equation (2,4), and minor losses 
(laminar) are included, the temperature rise will be 
increased slightly. The trend will be the same and the 
temperature rise at the higher RPM's will still be accep- 
table. The actual design of the closed-loop thermosyphon 
will require an optimization of the heat-transfer calcula- 
tions and the f r ict i cn- fact or calculations. This optimiza- 
tion may indicate that a larger, single, closed-thermosyphon 
loop for a conductor bar group is preferable over individual 
conductor loops. 

C- iDV4NTaG2S/DISADVABTAGZS OF LIQUID COOLING 

The main objectives cf cooling are to remove the heat 
created by the electrical losses and improve the operating 
efficiency, extend the lifetime of the insulation, and 

reduce the overall weight and size of the device. The 

feature cf liquid coding (vs. gas cooling) is the superior 
capacity of liquids, especially water, to remove heat with a 
small rise in temperature. The principle disadvantage of 
using liquids, or gaseous fluids other than air, is the 
design cf the piping, coupling, and secondary heat exchanger 
required to support the device cooled by the fluid. Since 
one of the main objectives is to reduce weight and size, the 
fluid chosen must be able to remove enough heat per unit 
mass to accomplish this goal within the device. 
Additionally, since the application under evaluation is 
marine combatants, the fluid should be non-toxic, and th<= 



68 




I I 



loe 




ON 



69 



Figure 2.7 Theraosyphon Analysis for Nusselt number. 



systaa should bs durable enough to withstand damage and 
maintain operational safely. The popular fluid to date ir. 
land-bassd installations has been hydrogen which is cl-arlv 
an unacceptable choice for the marine combatant due to its 
explosive nature and since the size reduction cbtain=’d is 
not sufficient to warrant system development. Of the 
remaining fluids, the obvious choice for its high heat 
capacity, availability, and safety, is water. The use of 
water has been tasted both in staters and, on a limited 



scale, in rotors. 

The principle problem to overcome in the use of water is 
in the rotary seals needed to direct the fluid into/cut of 
the machine. Additionally, when the system inside the 
device is analyzed, the ability of a design to withstand 
damage is questionable. The design would be critically 
dependent on the fluid within the rotor for operation, even 
in a casualty mode. The supply of cooling water to a stator 
would not be as critical, since any casualty could be 
repaired easily with shipboard damage-control equipment, and 
emergency cooling of a stator could b<= accomplished by 
external means due to the thermal conductivity of the stator 
housing. For the rotor, this leaves a choice; either not 
utilize the enhanced cooling and suffer the weight and size 
penalties, or develop a closed system for cooling the rotor. 
The closed system for the rotor is not necessarily isolated. 
A closed-lccp, thermesyphon could be placed in the rotor, 
transferring the heat from the rotor to the rotor shaft. 
This loop could have additional heat transfer to air passing 
through the core of the rotor. Further, the heat trans- 
ferred to the shaft could be removed to the ambient air or 
an external heat exchanger utilizing the non-coupling end of 
the mozor by the use of circulating water through the shaft, 
or by a rotating heat pips as seen in [Ref. 48]. 
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An addi'*‘ional method for cooling the rotor, which has 
bean studied by numerous authors, is the use ci heat rices 
or cicsed, tvo-phasa thermo syphons. In this case, the hsa-^ 
is transferred from the rotor to air circulated through the 
end-bell areas of the device. This method of rotor cooling 
greatly enhances the survivability of the device in a 
damage-control sense. Each heat pipe is independent of the 
others and damage to a few, unless all in the same group of 
adjacent conductor bars, would not jepardize the device; 
indeed, even if all were in the same group, limited opera- 
tion could continue. 
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III. TWO-P H &SS FLOiS COOLING 

The concept of utilizing the latent heat or vaporization 
to transfer heat is widely used in heat exchangers every- 
where^ frcm household and industrial boilers and space 
heaters to marine propulsion units. Air-conditioning equip- 
ment and heat-pump units also employ this idea to advanta- 
geously utilize the lew temperature of vaporization of 
f luorc-carbon compounds to remove the heat from the ambient 
and release this through condensation under pressure, thus 
transferring the heat in the desired direction at the cost 
of the energy to pump and pressurize the gas to the point of 
condensation. In each of these two-phase applications, the 
transport of the working fluid may require a pump, or it may 
be accomplished by the buoyant force that exists in an 
acceleration field, such as that 'due to gravity. This is 
the principle that is employed in the closed two-phase ther- 
mosyphon.- This device is similar to the classical "heat 
pips" and the only difference is that the "pumping" method 
in the heat pipe uses capillary action instead of accelera- 
tion forces to pump the fluid. Thus, a "heat pipe" used in 
a rotating reference utilizing the acceleration forces to 
transport the liquid from the condenser section to the 
evaporator section (in lieu of a wick) is, technically, a 
two-phase ther mosyphen. Figure 3.1 illustrates the opera- 
tion of the closed two-phase thermosyphon. The evaporator 
end must be in the direction of the acceleration field for 
the liquid to be transported from the condenser. Figure 3.2 
illustrates the operation of the heat pipe, which utilizes 
the capillary action of the wicking material to transport 
the liquid from the condenser to the evaporator. The 
remainder of the chapter shall use these definitions for a 
discussion cf both devices. 
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Acceleration 



F i 0 i d 




Figure 3.1 Typical Two-Phase Thermosyphon. 



Acceleration Field Not Required 




Figure 3.2 Typical Heat Pipe. 

A. HISTORICAL DEVEICFHEHT 

The development cf the heat pipe (which will apply 
equally to the development required for the two-phase ther- 
mosyphon) was begun by Gaugler [Ref. U9] in 1944 and by 
Trefethen [Ref. 50] in 1962. Trefethen was working on 

cooling methods for spacecraft (zero-g environment) for 
General Electric Ccmpany and discovered that capillary 
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pumping can be a very valuable area for further development. 
Working independently, Grover, en al. (Ref. 51] reinvented 
this concent ; he widely published rhe first applications of 
the device and gave it its name. Kraus, et al. (Ref, 52] 
have reported this historical development and have given 
sample calculations for the design of a heat pip*?. Chi 
[Ref. 53 3 also reviews the theoretical development of heat 
pipes, and provides the design procedures, including sample 
calculations, for their utilization. 

B. CONVENTIONAL, CLOSED TWO-PHASE THERMOSYPHON ANALYSIS 

As previously mentioned, the appealing factors associ- 
ated with the two- phase utilization include: high heat 

fluxes associated with phase changes, lower temperature 
gradien-^s associated with these processes, and reduced 
weight of the two-phase system over the single-phase liquid 
syst em. 

Research by Cohen and Bayley [Ref. 54] (referred to and 
discussed by Japikse [Ref. 44] ) found that the amount of 
liquid filling the system functionally affected to the heat 
transfer in both rotating and static tests. They found that 
heat transfer increased as the percent liquid in the system 
was increased to approximately 1.5% by volume, then 
decreased to an intermediate value, and finally increased to 
approximately the same value as in the 1.5% case. This was 
related to the following process: in the low filling situ- 

ation, condensate returns to the evaporator section and 
forms a thin film on the walls. It is in this film that the 
heat transfer occurs and the phase change takes place. In 
the case of the co mp letely- filled evaporator, a liquid pool 
region develops, and a condition of nucleate boiling exists. 
If the filling is insufficient, with regard to the heat flux 
being transported, the pool/film will not continue to the 
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end of nhs beared section and ”dry-out*’ will occur. This 
"dry-cut" will assume a drop-wise instead of film- wise chase 
change. The geometry of the regions in which the films are 
likely tc form, influences the results. 

Lee and Mital [ Eef . 55] conducted experiments with an 
electrically-heated, watsr-cooled thermosyphon using water 
and freon as the test fluids and varied the filling quan- 
tity, Lj^/L^, pressure (or heat flux (L^/L^ is the 

ratic cf the length of the heated section to the cooled 
section). The result of the filling quantity on heat 
transfer was the same as that reported by Cohen and Bayley; 
increasing heat transfer with filling to a point and then 
decreasing beyond that value to an intermediate value and 
increasing to the case of the completely-filled evaporator. 
The effect of decreasing l^/L^ was to increase the heat 
transfer within the range 0.8 to 2.0 and the advantage of 
larger condenser area was evident. The heat-tr ansf er coef- 
ficient was found to increase signif icancly with increasing 
mean pressure due tc at least three factors: 

1) since the mass flew for a given heat flux is nearly 
constant and density of the liquid increases with 
pressure, a lower pressure drop (and lower AT) is 
necessary for the same heat flux; 

2) for a larger varies much more rapidly with 

S O.U Sac 

T for the fluids considered, hence requiring smaller 
sat 

AT's at higher pressures; and 

3) for lower pressure drops, a more favorable force balance 
exists on the condensate film permitting a faster liquid 
return. 

Water was found to give heat fluxes superior to those of 
freon, for for the same aT, due to the larger values of 
latent heat of vaporation, and thermal conductivity. 



75 



Lae and Kital [ Bef . 55] also considered the ar.alY'ical 
problem cf pradicting the maximum heat-transfer rate for a 
laminar film, constant- w all-temperature ccndenssr and a 
constant-heat-flux evaporator. Neglecting the forces due to 
vapor pressure drop and momentum changes, and using a force 
balance on the falling film, balancing the effects of 
gravity and fluid shear, they related the mass flow rats to 
the heat flux. Using a local energy balance and an overall 
energy balance between the condensing section and the evapo- 
rator section, the following relations for the heat trans- 
ferred <q) and the saturation temperature (T ) were 

sat 

obtained : 

q = (3.1) 

an d 

g = {k(T -T )/(P.L, /L )}C/D (3.2) 

sc he 

where C = 1/8 - 1/2 ys + 3/8 y* - 1/2 y^ In y , 

D = y^ln y (1/2 In y2 -1/2) + 1/8 y^ +y2(i/2 In y -1/4) =1/8 

and y = 1 - 6/R. The simplifying assumptions (forces due 

tc vaper shear and momentum changes were neglected in the 
force balance) cause an error (of as much as a factor of 2) 
to occur, but qualitative behavior is correct as far as 
L^/L^ , working fluid are concerned. This develop- 

ment is also included in the survey by Japikse [Ref, 44], 
This heat transfer is graphically presented in the Lee and 
Mital paper in their Figure 12, shown here as Figure 3.3. 
This clearly shows that the quantity of heat that can be 
removed by a two phase water device is large, even at a low 
saturation temperature. The report also indicates that for 
water the increasing ratio of heated lengrh to cooled 
length, L+, decreases the maximum heat transfer. The heat 
transfer is also very sensitive to the operating pressure; 
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Figure 3.3 Compariscn of the Experiaental Results «ith the 
Analytical Predict-ion-f rom Lee & Wital paper. 



the heat -transfer coefficient increases appreciably with the 
mean operating pressure. The temperature drop increases 

with the pressure drop along the length of the thermosyphon 
tube (and the saturation pressure gradient). As the pres- 
sure increases, the vapor specific volume decreases, 
resulting in a decrease in the pressure drop; the mass flow 
rate of the vapor is essentially constant for transferring a 
given heat input rate. In the previous papers, the 

vertically-oriented thermosyphon was considered and the 
current model (for some of the orientations) requires the 
device to be oriented in the horizontal direction (perpendi- 
cular to the applied acceleration field) . Although the Lee 
and Mital paper was considering a vertically-oriented, two- 
phase thermosyphon, the general results are the same for th*= 
horizontal tube. 
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LiaiTS OF 0P2RATICH 



iihsn considsriag th? design of a two-phasa thsrn 
or a head pipSf liiitations to the heat transfer 
considered; four of thase limits are common to both 
(sonic, entrainment, toiling, and condensing). A 



os V phcr. , 
must. D = 
devices 
qualita- 



tive ccicparison of these heat-transf ar 
func'ticn of f or. f fs 

and is shown in Figure 3.U. 



limitations as a 
given in [Ref. h8] 




Figure 3.4 Operating Limits of Rotating Heat Pipes. 

The sonic limit and the boiling limit have been readily 
analyred. The entrainment and the condensation limits are 
not as well-known and little, if any, literature describing 
these phenomena exist. This is an area in which further 
research is needed. 

Sonic limii* vapor flow in a two-phase closed 

sysrem is limited to the sonic velocity at the operating 
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pressure (also known as the "choking" limit) . The sonic 
limi* is rspnesented hy: 



Q 

max 







A 

V 




(3 






) 



where is the sonic velocity of the vapor. This limiting 
velocity may be determined experimentally, computed by the 
relationship V ={(3r/3o) } , or apprcxima t«=d bv th'=' 

perfect-gas relationship {YET 

S 3. u 

2) Bnt r ein meet limit ; the interfacial shear between the 
liquid and the vapor will held the liquid (which is flowing 
in the opposite direction) back and starve the evaporator 
secticn cf liquid. This counter-current flow, when the 
relative velocity is large, causes the interface to become 
unstable which results in waves at the interface. As the 
vapor velocity continues to grow, droplets of liquid are 
formed at the liquid surface as the shear force exceeds the 
surface-tension force. The formation of these droplets and 
their subsecuent entrainment in the vapor stream causes the 
partial or total stoppage of the flow (dry-out). This’ 
phenomenon is generally governed by the Weber number (the 
ratio cf the inertial force to the liquid surface-tension 
force) . The Froude number is also used to characterize the 
phenomenon cf the drep-wise entrainment of the liquid in the 
vapor stream. In the application considered herein, the 
formation cf the waves is considered unlikely due to the 
extreme acceleration field present. Since the flow is 
counter-current and is in the presence of a high accelera- 
tion field the entrainment limitation thus will reduce to 
the "hold-up" limit. It will remain stratified until 
hold-up occurs, resulting in evaporator dry-out. 
Experimental determination of the exact correlation for the 
hold-up, or entrainment, limit is required for a horizontal 
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tube rc-atir.g around a parallel axis, thereby creating a. 
perpendicular acceleration field many times greater than 
cr^.vd.'ty* Th^ 9r. limit is G9n^rs.llv giv8n 

*^inax"^‘fg^v^v^E ' (3.4) 

where is when the Weber number equals unity, the Froude 

number equals unity, or is experimentally found. In another 
approach, Jaster and Kosky [Bef. 563 used the ratio of the 
axial shear force of the vapor flow on the liquid surface to 
the gravitational (accelert ional) body force upon the liquid 
(defined as ”F”) in order to arrive at an appropriate 
'*hold-up” criterion. They correlated experimental data as a 
function of this •'?” value, as the flow they were studying 
transiticned from stratified to annular. The result of 
their experiments was that the flow was stratified for ”F" 
values below five (5) and was f ully-annular above twenty- 
nine (29) . The value of "F" set equal to five (5) was 
established as the criterion for the transition from strati- 
fied to annular for the co-current flow case. The same type' 
of analogy was developed by Collier and Wallis [Ref. 57] who 
balanced the inertia force and the acceleration force to 
scale stratification effects where their criterion was given 
as j* (volumetric flux) equal to 0.25. Both of these refer- 
ences point out the need for further research in this area. 

3) Eoi lin q lisit; Again, as discussed in the first 
section of this chapter in the general discussion on two- 
phase thermosyphons, the boiling limitation in the current 
application is due to the creation of profuse nucleaticn at 
the evaporator section and a resultant vapor film on the 
evaporator surface which insulates the evaporator wall, 
resulting in dry-out and overheating. The boiling limit is 
given by the Zuber-Kutateladze prediction as discussed in 
[Ref. 48] and is given as: 
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( 3 . 5 ) 



Q =0. 13.^0 { o2P(p -p )a 

max V A.q 1 v 



vhare A 
haa t ing 



B 



is ihe area in tha evaporator section 



■f* T* ,'T 






4) Con den ser lia it ; The capacity of the condenser is 
dependent on many parameters, such as the geometry, nhe 
working fluid, the orientation of the vessel with respect to 
the acceleration field, and the operating conditions. The 
orientation of the vessel, for the model considered, is that 
cf a horizontal tube (perpendicular to the acceleration 
field) with internal condensation. The condensing liait is 
the capacity to condense the vapor on the inside surface of 
the condenser section, and this is shown in Figure 3.4. No 
analysis on the outer surface is presented in this thesis. 
Instead, the condenser limit is based on the capacity of the 
inside wall exposed to condensation, using an arbitrary, but 
reasonable, wall temperature. (Generally, the largest 
thermal resisrance occurs between the ambient and the 
condenser cuter surface. Therefore, heat transfer is 
limited by this resistance unless the outside surface area 
is substantially increased by the use of fins.) .The 
condensing limit is given by an equation from [Ref. 58] by 
Collier for a horizontal tube wirh internal condensation in 
a gravitational field as: 

Q =A ATF {(p (p -P )gh k 3) / 
max V c 1 ^1 fg f 



))■/ 



(3.6) 



For the rotating reference being considered, the "g” term in 
equation (3.6) is replaced by the centripetal acceleration 
cf The factor "F ” in this equation, allows for the 
fact that the rate of condensation on the statified layer of 
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liquid is r.«gl i cable , The value of *'? " depends on nh- 



angle d and is tabulated 


in fBef. 


58]. 


Its 


value ran 


ges 


from zero, when the tube 


is full 


(no 


surface 


remains 


for 


condensation to occur), no 


0.725 , 


when 


the tube 


empty. 


The 


angle $ is shown in Figure 


3.5. 












Figare 3.5 



Laminar Condensation within a Horizontal 
Tube (from Collier) . 



1 . Co m cut at io n of Heat -Tr ansfe r Limit s 

The calculation of these limits is an important step 
in the design process. 

The calculaticns required to analyze the configura- 
tion of the current model are not all well known. The 
following appear to be the most applicable to the model 
being evaluated: 

The sonic limit calculation is from basic principles as 

shown by equation (3.3) by allowing the vapor velocity to 
increase to the Mach value of unity. 

The boiling limit is given by the Zuber- Kutateladze 

prediction as discussed in [Bef. 48] and is given by equa- 
tion (3.5) 
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The condensing limit is given by an aquation frcni 

[Eef. 58} hy Collier in equation (3.6), 

The literature reveals no correlation yet for the entrain- 
ment limit in the heated, horizontal tube rotating about a 
parallel axis with counter- current flow. Since the majority 
of analyses previously conducted for this type of flow have 
shown that the entrainment limit is typically below th= 
sonic and above the condensing limit, the ccrreia-^ions 
previously mentioned have been plotted against the amount of 
fluid that is charged in the tube (percentage fill) for the 

three reasonably known values (sonic, boiling, and 

condensing limits) in Figure 3,7. Mso shown is the 

entrainment limit computed from the Jaster and Kosky 
[Ref- 56} correlation developed for co-current flow. 

As was previously noted, the condensing limit is th'= 
factor that is the most limiting in the design of the ther- 
mosyphon. The values of are increased to a maximum 

based on the length of the condenser, the percentage fill 
and the 1 T available between the constant temperature 
(assumed) wall and the saturation temperature of the vapor 
condensing on it. As already mentioned, the resistance to 
heat- transfer afforded by the condenser outer walls to the 
ambient may be g’lite significant. The ability to design an 
adequate fin arrangement and supply cooling air in suffi- 
cient quantities, while minimizing the windage losses, 
itself is a substantial problem. From the data obtained in 
this rough model, there are some self-evident points to be 
considered : 

The length of the condenser, from a heat-transfer point 

of view, determines the heat-removal potential, 

The amount of liquid fill can be optimized. 
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Tb-9 ten! per a t ure at which the walls of the condenser ar'^ 

maintained and the t emperat are at which the thersosyphcn is 
coerated are critical as this t'=m perature difference is the 
driving feres in the heat transfer. 

For the model being evaluated, and neglecting (as 
the limiting value) the entrainment limit, the closed, 
■^wc-phase t her mo sv phen, operated at a saturation tamperature 
at, or above, 100°C and with the condenser section held at a 
wall temperature of 50° C with rotation at a 0.381 m radius 
and 3,600 RFil can transfer the required 50 W of power with 
the condensing limit being the limiting value. Figure 3.7 

is based on a heat load of 50 W to be equally divided 

between condenser sections in both end-bell regions (i. e. , 
approximately 25 W each per condenser) . For the 3600 RPa 
model the percentage fill may increase to approximately 40%, 
and still maintain the required heat flux. The effect that 
the Er'1 has with regard to the temperature rise in the 
heated section, makes the choice of the closed-loop thermo- 

syphon poor for the motor application, but quite suitable 

for the generator. Further background on the theory of 
two-phase flow is available in publications by Wallis 

[Sef. 59] and by Collier & Wallis in [Ref. 573- The feasi- 
bility of the design has been shown by German and McLaughlin 
[Bef. 60] and by Sroll, st al. [Ref. 61] through experi- 
ments with heat-pipe cooling of electric motors in the 
scheme primarily recommended, although the difficulty of the 
liquid transport at lew rotation rates was noted. A diagram 
of the cooling scheme such as recommended and tested by 
Corman S McLaughlin and Groll, et al., is shown in Figure 
3.6 (reproduced from [Bef. 60] ) . 
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Figure 3.6 Typical Heat-Pipe Cooling Scheme 



Entraining Limit 
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Figure 3.7 Heat-Transfer Limits in the Closed 
Two-Phase Thermosyphon. 



cose LINING SOMMiRY 

As aiscussed in the paper by Greene, at al. (Ref. 62], 
the need for electric propulsion of !r.inimal cost, volume, 
and weight is well known and motivated by the inherent 
advantages of this type of propulsion. A significant 

constraint to be overcome in the development and manufacture 
of these systems is the rate of hear removal from these 
devices. In 1976, W. D. Morris [Ref. 63] mentioned than 
increasing the load per unit weight of these devices 
demanded that the inherent heat created by the internal 
losses be removed via improved cooling techniques. Morris 
also elaborated on the "complex arrangement" of the ducting 
and associated equipment required to accomplish this 

increased cooling. Further, he stated that "Ther<= is 
currently insufficient technical informarion available with 
which to make confident predictions of the effect of 
Coriolis and centripetal acceleration on flow resistance and 
heat transfer in the complex coolant circuits...". Since 
then, Morris and his co-workers have attempted to fill the 
literature gap in this area, and have uncovered a great many 
facts regarding this flow structure; they have also uncov- 
ered a great many questions for further research. 

As shown as early as 1955 in an article by Sir Claude 
Gibb [Ref. 64], the principle factor in the failure of two 
very expensive generating sets was the inability of the 
devices to dissipate their own heat. In the late 60 's, 
Feach [Ref. 65] wrote of the successes that the English 
Electric Company had with regard to liquid cooling of the 
large generator rotors being used in their 500 MW devices 
and predicted that the devices would reach 1,000 MW with 
this technology. This was followed in 1970 with a second 
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report tj Peach [Ref. 66] where he outlined the grcw-h o-^ 
D.S. generator sizes and predicted that the liqui d-cociing 
technology would be required in the near fu'^ure to supocrt 
the devices being planned in this country which exceaded 
1,000 as. The losses developed in the large-scale, high- 
power, superconductive devices, currently being developed 
for large-scale implementation, require extensive use of 
cooling at cyrogenic temperatures in order to function in 
the super-conductive stats in which they operate. A discus- 
sion of this application is presented by Schwartz and Foner 
[Ref. 67], in which a number of the schemes discussed herein 
are shown. 

The previous chapters have discussed the ability of 
external, forced-convective cooling to remove the 
internally- dev eloped heat loads in electrical, rotating 
machinery. The reduction in size, the increased efficiency, 
and the prolonged life cycle are all very important reasons 
to utilize advanced-cooling concepts. The use of turbulent 
flow is feasible at all PPM's; therefore, both moror and 
generator applications would benefit from forced-*convect ive 
liquid cooling utilizing turbulent flow. Laminar-flow heat 
transfer is feasible only at higher RPH's, which tends to 
limit the application of this method in regard to motors. 
The actual transition from laminar to turbulent flow is not 
well known and further investigation is required to deter- 
mine the transition point; the transition may occur at a 
flow and EPM acceptable for the motor application. However, 
the utilization of external, forced-convective cooling has 
some major difficulties associated with it; 

1) The source of the external cooling fluid must be 
easily accessible, dependable (in a casualty-control sense) , 
safe, volumetrically small, light, and should be reasonably 
inexpensive . 
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2) The fluid itself should be of high specific heat, 
light, inexpensive, and nore importantly, non-toxic and 
non-flammable. 

3) The method of transporting the fluid from the source, 

through an external heat exchanger, through the device, and 
back, should be extremely reliable and quire resistant to 
external damage. t his is an srea of s^ipus cpnsiderar ion 

rel iablility of conventional, r ot at in q sea 2^ is 
questi onab le in lea st. 

The alternative to the use of external, forced- 
convective cooling is the utilization of a closed device: 
(a) the closed, two-phase thermosyphon, (b) the closed, 
rota ting-loop thermo syphon , (c) the closed rot ating-loop , 

two-phase t her mosy phen, or (d) the heat pipe. These alter- 
natives remove the necessity of using radial rotating seals 
in the design of the cooling system. 

The closed, two-phase thermosyphon, as shown in Figure 
3.6, must have a sufficient forced-air circulation through 
the end-bell regions to remove the heat. ■ This will require 
optimizing the length of the condenser sections, the geom- 
etry of the fins, the number of fins and the direction of 
air flow within the device to maximize the heat transfer 
while minimizing the windage losses. The closed two-phase 
thermesyphon is feasible at the higher RPM's that the gener- 
ator operates at; additionally, this method could also be 
feasible in the motor application if a heat pipe (i.e., a 
thermosyphon whose walls are lined with a capillary wicking 
material) is used in lieu of the two-phase thermosyphon for 
the lower RPM's. 

The closed, rot ating-loop thermosyphon requires an addi- 
tional heat ‘exchange through the shaft of the device, 
possibly supplemented with forced-air convection through the 
hollow rotor region and by the exposed radial sections of 
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the loop. This shaft heat exchanger may be in the fcrm of a 
rotating he at- pipe, as shown in Figures 4,1 and 4,2, where 
the external exchange is to either ambient forced-air, or to 
an external, closed water spray. The shaft heat exchanger 
may also he cooled by forced convecrion using liquid 
cooling, with axial, packing-type seals, external to the 
device as shown in Figure 4.3. This closed-loop thermcsy- 
phon could use either single-phase or two- phase cperaticn. 

The closed two-phase th ermosyphon could also be extended 
along ehe shaft through a bearing plate arrangement and 
cooling air could he supplied external to the casing, 
simplifying the internal configuration, as shown in Figure 
4. 4. 

The ability of these devices to cool effectively has 
been proven and their usage in this application is theoreti- 
cally feasible. The exact correlations for the calculation 
of the performance values for these devices, especially in a 
rotating reference, will require further research. 

The analysis herein has shewn the ability of th® cooling 
schemes discussed to remove the .requisite heat load in order 
to extend the life of the insulation, increase the effi- 
ciency, decrease the weight and size of motors and genera- 
tors. This study graphically presents the applicable 

correlations for the determination of the heat transfer 
within a rotating, electric device. It alsc discusses the 
limitations of these correlations. It has been shown that 
closed, two-phase systems of an inherently higher reli- 
ability are feasible, and deserve further evaluation. 
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Figure 4.1 Closed-loop Ttiermosyphon to Shaft Rotating 
Heat Pipe (Air Exchanger) . 
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Figure 4.3 Closed-loop Theraosyphon with Shaft Forced-Convection Liquid 
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CONC IOSI ONS R.gCOMMENDaTION5 

a. CONCIOSIONS 

Eassd upon the information reviewed and analyzed in this 
thesis, the following conclusions are made: 

1) For ced- conv active liquid cooling is a practical 
means for improving the efficiency, reducing the size and 
weight, and extending the life of high-power motors and 
generators. The sain disadvantage of forced-convective 
liquid cooling is the requirement of radial rotating seals. 

2) Closed-loop liquid t her mosy phons are feasible at 

high BFM. They offer improved reliability over forced- 
convective liquid cooling which utilizes rotating seals, 
however, this method requires the use the shaft as a 

secondary heat exchanger. 

3) The closed, two-phase thermosyphon is feasible at 
the higher RPMs of the generators. This method requires 
forced-air convection through the end-bell regions or the 
extension of the devices along the shaft for external 
forced-air convection, both of which will add to the windage 
losses. 

4) For low RPM applications, such as motors, the use of 
a heat pipe will be required to overcome the inability of 
the low acceleration field to transport the cooling fluid. 
Both the closed, two-phase thermosyphon and the h^t-pipe 
require forced convection at the condenser ends +o remove 
the heat either internal, or external to the device. This 
may add to the windage losses. 
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3) Corrslat icr.s and theoretical analyses are net avail- 
able in the literature for: 

a> the counter-current entrainnssnt limit for t^o- 
phase ther nosy phens under rotation and 

b) the transition from laminar to turbulenr flow 
for rotating systems. 

E. BECCHMENDATIONS 

The following r ecemmend ations to continu=> this work are 
mad e : 

1) Testing of the forced-convective liquid cooling 
scheme should be accomplished using the DTNSHDC test rig to 
confirm the analysis and the validity of the available 
correlations, or to obtain proper correlations for the 
specific geometry and utilization being projected. 

2) An experimental program should be devised and 

conducted to determine the entrainment limitation in a hori- 
zontal, heated, channel rotated about a parallel axis with 
counter-current flow of liquid and vapor. 

3) An experimental program should be developed to 

determine the optimum condenser length and proper fin 
arrangment for the use of both heat pipe and closed two- 
phase thermosyphon systems for rotor cooling. Consideration 
of the structural integrity must be included for the 
condenser end-lengths. 

h) An experimental program should be developed for 

correlating the laminar- to- turbulenx transition during 
forced-convective liquid cooling in a heated, horizontal, 
rotating pipe. 

5) An experimental program should be developed to 

modify the inside geometry of the closed-loop thermosyphon 
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(ani axial, closed, two-phase theraiosyphor.) to optiniiz'^ the 
h-eat-transf =r and the friction losses within the device, 

6) A detailed analysis should be made of the shaft 
cooling potential for the rotating- loo p thermosyphon with 
regard to the heat-t ransfer enhancement cf the shaft to the 
fluid and the mechanical design of the shaft. 

7) An experimental program should be developed to 
analyze the closed, two- phase thermosyphon with radial 
condenser sections and the closed-loop two-phase thermosy- 
phon. These two devices may experience start-up balance 
pro blems . 

8) Various configurations should be analyzed for their 
potential for heat transfer to include a combination of the 
systems described herein, as depicted in Figures 4,1 through 

4.4 . 

9) The concept cf modifying the geometry of the rotor 
windings within these devices in order to minimize heat 
conduction path lengths to the casing, shaft, etc., and 
thereby the internal temperature gradient, should also be 
investigated. 
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APPENDIX A 

GOVERNING EQUATIONS FOR LAHI3AR CONVECTION IN ONIFOEMLY 
SEATED HORIZONTAL PIPES AT LOW RAYLEIGH NDHBERS 

The followiag development, by Norton [Ref. 11], ie 
presented herein for completeness. The steady laminar 
motion of a fluid in a horizontal circular pipe of radius a, 
the walls of which are heated uniformly so that a constant 
temperature gradient x is maintained in the direction of 
the axis. The flow will be referred to in cylindrioal coor- 
dinates (r,p , 2 ) with (J) measured from xhe upward vertical 
and the z-axis along the axis of the pipe; the velocity 
components are denoted by u,v, and w. The effects of dissi- 
pation and of the pressure term in the energy equation will 
be neglected and that variations in the density due to the 
temperature differences are so small that they only effect 
the buoyancy term. Thermal conductivity and kinematic 

t 

viscosity are assumed to be constant (which will introduce 
quantitative errors into the solution, but should not change 
its general character) . 

The equation of continuity is then: 

8ru/3r + 9 v/ 9 <{) + 9rw/9z =0, (A.1) 

and the energy equation is; 

u(9T/9r) v/r(9T/9(J)) + w (9 T/9 z) = k V2T, (A. 2) 

where V2=92/9r2 +1/r(9/9r) + 1 /r 2 (3 2 / 34 ) 2 ) + 32/922 and T is 

the local fluid temperature. The momentum equations can be 
written in the form of equations (A. 3) through (A. 5). 
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u ( 3u/oi) +v/r ( 3u/S<p) +w ( 3U/3Z) “v2/r = -1/p(3p/3r) 



v( V2u-u/r2-2/r2 (3v/3(i>)- g (T -T) cos.,*; , 



w 



l ^ • 3 ) 



u ( 3v/3r) +v/r ( 3v/ 34>) +w ( 3v/ 3z) =uv/r = - 1/pr ( 3p/ 34>) 




u ( 3w/3r) +v/r ( 3w/3(t)) +w ( 3w/3z) 



-1/P (9p/3z) +v ySw, 



(A. 5) 



In equations (A. 3) and (A. 4) , the buoyancy force has been 
calculated relative to the fluid at the same level adjacent 
to the pipe wall and the remaining distribution of force has 
been absorbed into the pressure, p. 

For steady convection sufficiently far from the pipe 
opening to avoid inlet-length effects, the temperature 
throughou-^ the pipe increases uniformly with the distance 
along the pipe axis. Hence, the dis"tribution of buoyancy 

force in sections of the pipe is independent of z; as the 
secondary flow is caused by the buoyancy forces, the flow 
field must also be independent of z. It follows that there 

should be a similarity solution with u,v,w, and T, -T as 

w 

functions of r and 41 only, and with p=Yz + P(r, 4> ) (where P 
contains the terms absorbed into p) . 

The continuity equation reduces to the form: 



hence, a dimensionless Stohes stream function ^ can be 
introduced in such a way that: 



3 (ru) /3r+ av/d<p = 0; 



ru/v =d^ / 3(f)r 



V/ V = - dtp/ dr 



99 



Ths main aquations (A. 1) through (A. 5) may ba racuced to 
non -dinisnsicnal form by tha transformations r=aH, w=(v/a.}!^, 
and T -T=TaPr , whers ?r is the Prandtl number (v/a.) • If 
the pressure is eliminated between equations (A. 3) and 
(A. 4), the resulting momentum equations are: 

7t + V? OVSHO/at) -3'i>/36 0/9?.) } v2;i; = 



5a {(90/3?) sine}) + 1/R (30/9 (t>) cos <j)} , (A. 6) 



V2W+ 1/E{(9>K3H) 0/9*)- 



(3*/3*) (3/3R)}W +4Re =0 , 



(A. 7) 



and where, 

Vf = 3 2/3R2 + 1/R0/3R) + 1/R2( d ^/d<p 2). 

The Rayleigh number, Ra, is given by g6Ta*/(av ) and the 
Reynolds number is the normal Reynolds number. Re, for 
Foiseuille pipe flow based on the pipe diameter and the mean 

velocity across a pipe section is (a/v) ( )- 

The solution of the energy equation depends on the temp- 
erature on the boundary condition at the wall of the pipe. 
Pipes will usually have re asonably-thick walls of material 
with thermal conductivity much higher than that of th<=- 

fluid, so there will be little variation in wall tempera- 

ture around the pipe circumference. This will be specially 
so for slow rates of heating when asymmetries of the flow 
will have small amplitude. Hence, write T =T0+ z, where TO 
is the wall temperature in the section containing the 
origin. Although there is uniform heat transfer per unit 
length of pipe, the local heat transfer will be slightly 
greater near the bottom of the pipe than near the top. With 
this assumption, equation (A. 2) reduces to the non- 
dimensional form of equation (A. 8). 
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& 




V20 + (a/H) [(3’V9R) (3 0/9<1>) -(9iV9<t>) (96/0R)}+W=O 



(A. 8) 



The houndary conditions are: 

\i,7fV, and 3 are zero on r=a; 
u,v,w, and 0 are finite on r=0, 

(A. 9) 

Although solution of equations (A. 6) through (A. 8), satis- 
fying conditions of (A. 9) is considerably difficult, succes- 
sive approximations to the solution can be obtained by 
expanding W, and 0 as power series in the Rayleigh 

number, provided that this is numerically small. Supposing 
tha t : 

i/;=Ratl^j+Ra2i{; +. . ., 

W=W g + SaW j + Ra 2 '*’ • • • r 

0=0o + Ra0i+Ra^02'*’.«« . 

(A. 10) 



The leading term of \p must vanish because there is no 

circulation when A is zero, but 9 q is not zero since the 

difference in temperature between a fluid element and the 

neighboring wall is proportional to t 0 • When relations 

(A. 10) are substituted in equations (A. 6), (A.7) , and (A. 8) , 

three sets of equations for the funtions \p , vj ., and 9. 

i 1 1 

are obtained by equating coefficients of powers of Ra. 

From equation (A.7), the basic equation is: 



V2«o + 4Ee = 0, 
which has a solution: 



(A. 11) 



Re(1-r2) 



(A. 12) 



which satisfies the conditions M g=0 at R=1, Wg finite at 
E=0. This is ordinary Poiseulle flow under a pressure 
gradient -UpN^Re/a^ in an unheated pipe. 
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The first of the equations derived from equation (A, 3) 
for the temperature distribution is: 




subject to 9g=0 at 5 = 1, 9g finite at R = 0, and the solution 

is axially symmetrical as no account has been taken of 
gravity at this stage of the approximation. This is satis- 
fied by: 

6 o = (1/16) Ee( 1-R2) ( 3-R2) , (A. 14) 

which is the customary solution for forced convection. 

The first-order approximation for ^ satisfies the equa- 
tion: 



7*ii^i=Oeo/a3) sin (i) (A. 15) 

obtained as the coefficient of Ea from equation A. 6 , and the 
boundary conditions 34 ; i/ 3R , i /34>=0 on R=1, and R~ ^ (94' 1 /3 j>) , 
34 )i /3 r remain finite at R=0. If 4'' is assumed to have the 
form 'j^j(R)sin4> , the dependence of equation A. 15 on 4) is 
eliminated, and 4 ;j(R) is easily found: 

=- (1/4608) ReR ( 1-E2) 2{ 1 0-R2) sin 4) . (A. 16) 

Similarly, the first-order approximation to W^satisfies the 
equat ion : 

V2W = (1/R) (34' i/3 4)) (3Wg/3R) , (A. 17) 

and the boundary conditions, M^=0 at R=1 and Mj finite at 
R=0. The dependence of equation A. 17 on 4> is eliminated by 
taking =W ^ (R) cos 4 ) , whence 
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W ^ = - (1/18U32 0) Re2R (1-R2) (49-51R2+19R^-R6) ccs<^>. 



(?-. 18 ) 



The fir.9t-ordar approx iia'icn ro 0 satisfies 



tion 



V^e = (o/K) (dip/d4>) (3 9 /3R) -H 



1 



(A. 19) 





ons, 9, =0 at 3=1 and 9 finite at 



3=0. Hence, 

0 1 =- (1/22118400) 89 2E(1-R2) { ( 38 1 + 1 325a ) - ( 354 +1 6 7 5a ) E2 ♦ 



This ccinpletes the approximate solution to the first 
order in Ra for convection in a horizontal pipe, which is 
heated so that there is a constanr temperature gradient 
along the walls and uniform temperature around the girth. 
It may he noted from sgua tions (A. 10), (A. 12), (A. 14), 

(A. 16), (A. 18), and (A. 20) that the full convection solution 
depends essentially dr. the product RaRe. This is clear if 
it is recalled that Ra is proportional to the increase in 
wall temperature along the pipe and Re is proportional to a 
characteristic velocity of the flow along the pipe. An 
increase in Ra means that the fluid will be carried through 
a larger temperature differences, and, hence, there will be 
increased buoyancy forces; but the same effect can be 
produced by an increase in He, and hence in the velocity of 
the main flow. 

The two most interesting consequences of such a solution 
are the modification of the heat transfer due to secondary 
flow, and the effect of heating on the tangential stress at 
the wall (or on the rate of flow for a given pressure 
gradient). However, as both W ^ and 9^ vary with cos <J> , 



+ (146+S2 50 ) R^-(29 + 200a) R*+(1-2a) R8}coS({) 



(A. 20) 
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neiTher the flux ncr the 
section of the pip's will 
apor csina ti cn a 1 th cuoh than 
To find ■‘■he changes ir- 
is necessary to proceed to 
and the revelant equations 



heat transfer across the yhole 
ch2.ncT 3d bv if irsh'~cr-?zr 

s will be local variations, 
tctal flux and heat transfer, it 
the second-order approximations, 
for these are; 



v^ :.s = (V-R) {3( ,Hj ) /d(R ,0 )} 

♦(38j/9 R) sin,j)+ (1/R) ( 3 e/ 3 c{)) coscj) , 

=(1/R) {S('« )/3 (R, <{))}+(1/R) (3(W , J/3(R,4))} , 

12 11 U ^ 

V^g =(a/F.) {3(9^,<|;^)/3 (E,({)) } + (a/R) {3(9^,ii^^/3 (R,$) 

(A. 21) 



The only difficulty 
tediora and the 
presented in [Ref. 



in solving 
r emaining 

11 ]. 



these equations is numerical 
second-order solutions are 
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APP^DIX 3 

SINGLl-EHASS ANALISIS (COMPOTER PROGRAH AND RESULTS) 



The following program was written to give the results of 
the various correlations for the heat transfer for th“ 
forced/ single-phase convection in she rotating reference as 
disscussed in Chapter II. The program is so written that 
modification to include other correlations is quite easy. 
In was written in HP-Basic for the HP-9826 computer. 
Sub-programs for the graphics are not shown. Sample runs 
are included immediately following the listing. 

-Pages 105- 109 Program Listing 

Therraophysical Properties (Common to 
Appendices B through D) . 

Sample Calculation 



-Pages 110-111 
-Pace 112 
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I 




Uo'J ' 
■0-40 
1050 
1050 
*070 
- 08 C 
1090 
; 100 
; ‘ 10 
1 120 

• *30 
1 1^0 
1 ’50 
1150 
’ ‘70 
1180 
1 ’90 
1200 
’210 
1220 
1230 
12^0 
^250 
1260 
’270 
1280 
’ 290 
1300 

• 310 
-320 
! 330 
1340 
‘350 
•360 
•370 
•380 
•330 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
-500 
1510 
1520 
1530 
1540 
1550 
•560 
’570 
1580 



“ ■fr 4 '’^h. ^ r^'. R 

CR2A:^l': janiiary 13. ’984 



COM /CiTun/ Xmin.Ynin.Sfx.Sfy 

printer is 1 



jSING “2X . "•'D'?-^* au-I t values: 

n=y 7ror_o ♦ ^ / 



D-4.76 


3E-? t 


Tuoe a^anerer 


( m ) 






L- .823 


0 r 


^uoe length 


(m) 






R- .402 


'i 1 


Radius of rotor 


( m ) 






PRINT 


USIfJG '•'IX, 


,*“'TuDe diameter 


= 


••••.Z.4DE.’^^ 


( m ) : D 


pRTNT 


USING ‘'AX 


.••“Tube length 


= 


••••.DD.D.^’^’ 


< n ) : 


PRINT 


USING ••4X. 


• I • t ^ _ 

. rvacius or rotor = 


••••.Z.3D,"‘^ 


(h) :R 


BEEP 












INPUT 


'•OK TO ACCEPT DEFAULT VALUES 


( 1-Y.O-N)?”, 


.Id 


BEEP 












INPUT 


"LIKE A HARD COPY ('=Y.C= 


N)?- 


. Ihc 





IF Ibc=1 TAEN 
PRINTER IS 701 

?rp1ter is 1 

END 

1“ Id-1 THEN 1300 
BEEP 

INPUT ••ENTER TUBE DIAME*"ER (.'n)".D 
BEE^ 

INPUT ••ENTER ^UBE LENGTH 
BEE'=’ 

INPUT ••ENTER RADIUS OF ROTOR 

^RIN"'^ USING •• 1 OX . ••••♦*■*' Geometric V5.riables 

^^'INT USING •• • 4X . ••••TaDe diameter <D) - “•• . Z . 3DE . •“• <n) :D 

PRIf'lT USING •• 1 4X . ••••Tiibe lenqtn <L) - •*" . DD . DD . “'• (m) :L 

^RINT USING •••4X,-’^L/D ' - **•• . 5D . D‘^ :L/D 

PRINT USING ••14X. ••••Rotor radius (R) - ’••• . Z . 3D . ’•’• <m) :R 

^PINT 

“RINT USING •• • OX . ••••ComDLLtec v/a.u.es are: 

PRINTER IS 1 
BEEP 

^RINT USING ••ZX. ■•••SELECT OP'^ION: 

PRINT USING •'4X.****- Sinaie point 

PRINT USING ••4X.*’”2 Nu Cersus RPN 

PRINT USING ’•4X/^^’3 Nu. versus Re 

PRINT USING ••4X/^*^4 Nu versus Prod 

INPUT Op 
IF 0p>3 THEN 

PRINT USING **4X.^*^^SELEC': OPTION: 

PRINT USING ••4X/^^’2 VARY RP^ 

°RINT USING ••4X/*‘^3 VARY Re 

INPUT Opx 

IF Odx-2 ^HEN 0p=2 
IF 0px= 3 THEN 0p=3 
END IF 

IF 0p>1 THEN 
BEEP 
J ^ = 0 

=>RINTER IS 1 

°RINT USING -4X. ••••SELECT OPTION: 

PRINT USING ••4X.“" 1 =Di ttus-Boel ter . Z-NaJiayama. 3-Nakayama d 



Fuz 1 ot a' 
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1590 


PRINT USING 


mar , 


8=Na>: avana” I 


1600 


ir^PUT Ian 


1 C 1 


r-vt-'-rk 


i oZu 


imroT '!_I\E Ti 


1630 


IF Ok'plot-1 


1640 


CALL Plot(Jj 


1650 


BEEP 


1560 


INPUT ’’SOLID^ 


1670 


END I'F 


1 6S0 


END IF 


1690 


Ns tep= 1 00 


’700 


I" Lip = 2 THEN 


1710 


BEEP 


1 720 


INPUT ’’ENTER 1 


1730 


Rpm=Rpm 1 


I 740 


lLSE 


1750 


BEEP 


1 760 


INPUT ’’ENTER i 


1770 


END IF 


1780 


BEEP 


1790 


INPUT ’’ENTER 1 


1800 


IF Ihc'1 THEN 


1310 


IF Ihc-0 THEN 


1820 


PRINT 


1330 


^RINT USING ” 


1840 


PRINT 


1850 


PRINT USING ” 


I860. 


IF Io=l THEN 


1370 


BEEP 


1380 


INPUT ’’ENTER 


1890 


BEEP 


1900 


^RINT USING ” 


1910 


INPUT ’’ENTER 


1920 


^RINT USING ” 


1930 


PRINT USING ” 


M40 


END IF 


1950 


IF Io=2 THEN 


1960 


BEEP 


1970 


INPUT ’’ENTER 


1980 


PRINT USING ” 


’990 


BEEP 


2000 


INPUT ’’ENTER 


2010 


PRINT USING ” 


2020 


PRINT USING ” 


2030 


END IF 


2040 


IF Io-3 THEN 


2050 


IF Jt>0 then 


2060 


3EEP“ 


2070 


INPUT ’’ENTER 


2080 


PRINT USING ” 


2090 


BEEP 


2100 


IF 0d<3 then 


2110 


INPUT ’’ENTER 


2120 


PRINT USING ” 


2130 


ELSE 


2140 


INPUT ’’ENTER 


2150 


Mf = Mf I 


2160 


END IF 


2170 


BEEP 


2:80 


INPUT ’’ENTER 



U ruui 
THEN 



^ = />iooQs d rloms. 5 = S teohenson . 6 = Woods*ifior r is- ^ 

V i = [.u = N> ^ ‘,Ukpiot 
. X , Y . Type , Cx > 



= L. am 



lOX.””*** Operating ^^ariabies * 
10X,*”*Input variables are: 

COOLANT MASS FLGU RATE (kg/s)”, 

MX . ””Coo lant mass fioui rate = 
INLET^^AND OUTLET TEMPS (C)”.^i. 
MX. ’"’Coolant inlet remp 
MX, ’’’’Coolant outlet temp 



IMPS (C>” 



””,Z.3DE. 

To 

””.DD.DD. 

””.DD.DD, 



MD.DD, 



’’.DD.DD, 



( k g/s ) ' 

(C) 

(O* 



;Mf 



:Ti 

':To 



(N)’ 



(O’ 



•’’.DD.DD,”” (O’ 



MD.DD.”” (W)’ 



:Ti 

;To 



■;0 



””.ZMDE,”” (kg/s) :Mf 
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2^30 

«i2c 

2230 

22A0 

2250 

2260 

2270 

2280 

2300 
2310 
2320 
2330 
23^0 
2350 
2360 
2370 
2380 
2390 
2400 
2^10 
2420 
2430 
” : To 
2440 
2450 
2460 
2470 
Ta 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620! 

2630! 

2640 

2650 

2660 

2670 

2580? 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 



PRINT USING *' 1 4 X . ”*’Coo i.ari t inlet temp = “”,DL.DD.”” <C> 

= " ^20 
END 

I a- V I o -r T 1 > ^ . 3 
Up=h NUpu ( Ta) 

Ma=FNriuu ( Ta) 

K-FNKu(Ta) 

Pr=FNPru(Ta) 

Rho^FNRhou ( Ta) 

Nuv=Mo/Rho 
IF THEN 

Q=Mf-*CD^(To-*l’i) 

PRINT USING ” 14 X.””Heat load <Q) = ””. 40 . D,”” (W) ;Q 

END IF 

IF Io -2 THEN 
Mf = 0 / (Cp*( To“T i ) ) 

PRINT USING ’* 1 4 X , ’*”Coo lant mass fiou rate = "** . Z . 4 DE , **” <kg/s)‘ ;Mf 

END IF 

IF Io =3 THEN 
Toc = Ti-*- 0 / <Nf*CD> 

IF HBS(To-Toc )> .01 THEN 
To= < To-^Toc ) *■ . 5 
GOTO 2220 
END IF 

IF Jj -0 THEN PRINT USING ” ^ 4 X . '*”Cool ant oiitiet temp = ’‘"» 3 D.DD,”” (C) 



END IF 

IF Jr=0 THEN 
PRINT 

PRINT USING ” 1 0X , "*T 1 LLid properties evaiaated at ‘*”.DD.DD, ’*'* (C) are: 



=^RIMT USING ” 14 X.’**’SDec:f ic neat (Cp> 

PRINT USING ’* 14 X. ••“Viscosity (Nu> 

PRINT USI^^G ■• ’ 4 X , ’•••Tne rmal cond (V) ~ 

PRINT USING •* 14 X .••••Prandti niimber (Pr) 

PRINT USING “ 14 X. ••••Density (Rho) 

PRINT USING •M 4 X. ••••Kinematic vis <Nuv) 

Be ta=^NBe ta( Ta) 

PRINT USING “MX.^^^Coef t.her exp (Beta) 
PRINT 

PRINT USING “ 1 OX . *•••♦♦♦ CaicLilat ions 

PRINT 

PRINT USING •• 1 OX . ••••P re 1 im 1 nary calculations 
END IF 

Re= 4 *Mf/(PI*D*Mu) 



•••• (J/kg.K) :Cp 

••••.Z. 4 DE.‘^^‘ (N.s/m' 2 ) :Mu 

■••• .Z. 4 D,““ <W/m.K) ;!< 

•••\Z. 3 D“:Pr 

•••• , 4 D.D.“^* (kg/m 3 > :Rho 

•••• . Z . 4DE , •••' ( rr 2/s ) ; Nuv 

••••.Z. 4 DE/*^’ < 1 /K) :Beta 



Friction Factors for stationary reference 
IF Re<2*1 .E+4 THEN 2670 
F=. 184/Re^ ,2 !IdD Eqn 8.21 
GOTO 2530 

F=.3l6/Re^.25 !IdD Eqn 8,20 

Qpp=Q/(PI*D*L) 

0p=Opp*PI*D 

Vm=4-rif/(Rho*PI*D'2) 

Dps*F*Rho*Vm"2-*L/(2^D> 

0mega=Rpm*2*PI/60 
Ro-Om/ (Omega^D) 

Jay =0mega*D'2/Nuv ?From STEPH-Rotat lonal Re 
Kil'*FNKui( Ta) 

Dt t* ( To-T 1 ) /L 
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27S0 Gr = ^^ ( Jnega 2 > *Be ta-* < ( D/2 J ^)*Dtt/Nn.v 2 
2790 Ra=Gr*Pr 

2c 00 Ganria= ( Re ' ( 22/ I 3 > > "" ^ ( Gr* < ^ r . 6 ) ) ( - I 2/ I 3 ) ) 

r> / i r- C\ 

C.O t \J V OOjKI.ia. w . ^ 

2^i52U ! 

2830! Laninar Value 
2840 Nular^-48/11 

2850! 

2860! Turbulent Value from D 1 1 tus-Boe 1 ter 

2870 Nudb».022^Re^.8*Pr* .4 

2330! 

2890! Steonenson Correxation 

2900 Nus= ( ?r " . 4/ . 367 ) ♦ . 007 ; *Re ' . 38* Jay . 023 ! Pi r , Turb“EQN2 . 23-cor rected for uat 
er 

2310! 

2920! Nakayana Correlation for turouient conditions 

2930 Nunak*(Re''.8)*<Pr^.4)*(Tg' (l/30))*(l + .014/(Tg)* (l/6>)*.033 ! Eqn2.9 <Tur 

blent) 

2940 Ty«Nunak / ( Re" . ‘ . 4 ) 

2950! 

2960! Nakayama Correlation for laminar conditions 
2970 J2-J-2 

2980 $q=2/n *n-^S0R( l+77/^-( l/Pr'2) ) ) 

2990 Cf = 1- .486*( (3*Sq-l >" .4/(Sq*(Sq*Pr*S0R(5)+2) ) >*J2/(Ra-Re) ' .6 
3000 Nunak 1=48/1 1-. 1 9 1 /$q*( 3^Sq- 1 ) ' . 2* < Ra*Re ) ^ . 2* 1 /< 1 +< I /10-Sq*Pr > ) 

3010! 

3020! Nakayama/Fuz loka Correlation for Radial Pioes 
3030 Nunf=( .014/.023)*(Re/Ro'-2.5)' .12^*Mudb ITURB-EQN 2.21 
3040 Mus = hu ITemp Value Mus’^Mu 
3050! 

3060! Sieaer-Tate Correlation for Turbulent 
3070 Nust=.027»Re" .8*Pr" . 3333* ( Mu/M us ) .14 

3080 Prod = Ra*Re-^Pr 
3090! 

3100! Woods-Mor r IS Correlation for Laminar 
3110 NuuP=.262*Prod' .^73*48/r. ! Ean2.22-WM3 
3»20 J1=Jay/8 ‘ 

3130! 

3’40! Wooas-Morris Correlation for Radial ^ipes 
3150 Nuum2= . C 1 5*Re' . 78* J1 " . 25 !from um2 
3160 J; = 0 then 

3170 PRINT USING " 1 4X . ’"’Reynolds number <Re) 

3180 PRINT USING ’’1 4X . ”’’Fr ict i on factor (stat) 

3^90 PRINT USING ”i4X.””Heat flux (Qpo) 

3200 PRINT USING ”14X,*’”Mean fluid vei (Vm) 

2210 PRINT USING •’14X,””RPM 

3220 END IF 

3230 IF lhc=0 THEN PRINTER IS 1 
3240 PRINT 

3250 IF Jj = 0 THEN PRINT USING ” 1 OX . ‘"’Resul ts : 

3260 IF 0kpiot=0 OR (Okpiot=l AND Nstep<11) THEN 

3270 PRINT USING ’MOX,”” Nudb Nunak Nunf Nuum Nus Nuuim2 
3280 PRINT USING ’M 3X . 6< 3D . DD . 2X )” ;Nudb . Nunak . Nunf . Nuum . Nus . Nuum2 
3290 PRINT USING ’M4X/"* Ra Ro Ra*Re*Pr Mf Nak 1 ” 

3300 PRINT USING "1 4X . 4( Z . 2DE . 2X) . 3D . DD , DD . DDD . Z . DD” : Ra. Ro . Prod . Mf » Nunak 1 

3310 PRINT USING ’M4X.”” RaRe Tg Ty 

3320 PRINT USING ”1 3X . D . 2DE .2X . Z ."3DE . 2X . Z . DDD” : Ra*Re . Tg . Ty 
3330 END IF 
3340 IF 0kplot=1 THEN 
3350 IF Ian=1 THEN Y=Nudb 
3360 IF Ian=2 THEN Y*Nunak 



= ••”.Z.3DE’’:Re 
- . Z . 3DE” : F 

= ”” ,*Z.*3DE,’’” (N/m'2> iQpp 

= ’"’.Z.2DE.”’’ (m/s) :Vm 

= ””,^D.DD/’:Rpm 
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3370 

3380 

3390 

?6Qn 

3410 

3420 

3430 

3440 

3450 

34S0 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 



IF Ian=3 THEM Y=Nunf 
IF Ian=4 THEN Y = [>iaiiMn 
IF Ian=5 THEN Y=Nus 

T- V = ,\,i u.irr, 7^ 

IF Ian=7 THEN Y=Nalam 
IF Ian=8 THEN Y = Niir.ak 1 
IF Op'2 and 0px<>2 iHEM X=Rpm 
IF 0p=3 and 0dx<>3 THEN X=M-r 
IF 0px=2 OR 0px=3 THEN X=Prod 
Jj=Jt+1 

CALL Plot(Jj.Jc.Jd.X.Y.Type.Cx) 

END IF 

IF 0p=2 THEN 
Rpm=Rpm*1 0~ (Cx/Nstep) 

IF RpnXRpmh THEN 3590 

GOTO 2730 

END IF 

IF 0p=3 THEM 

NT =Hf* 1 0' (Cx/Nstep) 

IF Mf>Nfh THEN 3590 
GOTO 2040 
END IF 
BEE? 

INPUT "ANOTHER RUN ( 1 = Y . 0=N ) P" . I r 
PRINT "PU" 

IF Ir=I THEN 1350 

INPUT "WANT TO L ABEL"? ( 1 =Y . 0=N ) " . 1 1 

IF Ii=1 THEN CALL Labei 

END 



110 



"JOO DEP :5am> 

• Ci 0 Dir* k<6> 

•om;i r>arci -7 ^ • I'T-'Sd-J . -It “ I . ' 7 '-c/.p 

un.M M.io/*ii7o2t2C.373Go/u.^.c7.t* 

!040 READ K(') 

’ 050 T=<Tsteajn+273.15>/647.3 

) 060 Smn=0 

1070 FOR M=0 TO 4 

1 080 Sum=SujTi+K(N)*( I-T) '(N’-l ) 

1090 NEXT N 

1 ■ 00 Br = Sun/< T*(iT<(5)-<'-T)+<(g>«<’-T)-2))-(1-7)/(K(7)*(1-')2+<(8)> 
1110 F r *EXP( Br ) 

'120 P=22120000*Pr 

1’30 RETURN P 
1140 FHEND 
' '50 DEF FNHfar^) 

I'SO Hf<3 = 2477200-2450’<T-10) 

1170 RETURN Hfc 
1180 FNEND 
1190 DEF FNMuu(T) 

:200 9=247. 3/(T+133. 15) 

1210 «u=2.4E-5-10''fi 

1220 RETURN Ma 
1230 FNEND 
1240 DEF PNVvst(Tt) 

1250 P=FNPvst<Tt) 

1260 T=Tt+273.l5 

1270 X=1500/T 

1280 F1-I/( !+T*l .E-4) 

1290 F2-(1-EXP<-X) )-2.5=EXP(X)/X' .5 

• 300 B= .0015*F-- .000942*F2- .0004882*X 

1310 K=2-P/(461 .52’T) 

•320 V=( 1 +( 1 +2'B*K) ‘ .5>/K 

1330 RETURN V 
• ?40 FNEND 
1350 DEF FNCDua) 

•360 Cdu= 4.21 120358-T-(2.2632oE-3-T-(4,4236lE-5+2.7l428E-?-T) ) 

1370 RETURN CdwIOOO 

1380 FNEND 

1390 DEF FNRhoui(T) 

1400 Ro = 999.52946 + T*< . 01 269-T-(5 . 4825 1 3E-3-T- i .234147E-5)) 

1410 RETURN Ro 
1 420 FNEND 
1430 DEF FNPrw<T) 

1440 P nii'FNCoiii ( T)=FNMuw<T) /FNKu ( T ) 

1450 RETURN Pru 
1460 TNEND 
1470 DEF FNHw(T) 

1480 X-(T+273. 15)/273. 15 

1490 Ku.--.92247+X-<2.8395-X-( 1 .8007-X*( .52577- .07344-X) ) ) 

1500 RETURN Ku 

*510 FNEND 

1520 DEF FNTanh(X) 

1530 P=EXP<X) 

1540 0-EXP<-X) 

1550 Tann=<P-0>/(P*Q) 

1560 RETURN Tanh 
1570 TNEND 
1580 DEF FNHf(T) 



1 1 1 



’S90 ^f=T-(^.2038^9-T-^5.38132E-^-T-^.551603^7 

1600 RETURN Hf*1000 

^610 -NEND 

1^20 DE" FNTx/sr'<f^' 

' o3G I i 1 0 
1640 TU1C 
i 650 Ta= < Ta+Tl )♦ .5 
"660 Pc*FNPvs t < Ta) 

1670 IF ^^BS( (P-Pc>/P)> .001 THEN 

1680 IF Pc<P THEN Tl=Ta 

1690 IF THEN Ta=Ta 

1700 GOTO 1650 

1710 END 

1720 RETURN Ta 

1 730 FNEMD 

1740 DEF FNBeta(T) 

1750 Rop=FNRhou ( T+ . n 
1760 Rom=FNRhow<T-.i ) 

W70 Beta='-2/(Rop+Rom)-^(Rop-Rom)/ .2 
1780 RETURN Beta 
• 790 FNEND 
1800 DEF FNfilpha(T) 

1310 Alpha=^NKu< T ) /(F NRhouiC T )*FNCpu( T ) ) 

1820 RETURN Alpha 

1330 "NEND 

1840 DEF FNSigtr.a<T) 

1350 Tt-547.3*-T 

I860 A*.00WTt^2*( . 1 1 60936807/( 1 + . S3*Tt ) ) 

1870 B- .001 121 404688-5. 752805 18E-6-Tt 
^380 C=i .28627465E-8*Ti 2- 1 . 1 497 1 929E- 11 -T t '3 
1890 S iana*A'^B'*'C 
*300 RETURN Signa 
, 1910 FNEND 



E - 6 ) ) 
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f;^C!netric /ar:aijlef 
ciiaiReter (D) = 

L/ u 

Rotor radiLLs <R> 
Conotited valaes are: 



4.7S3E-03 (n> 

no / . > 

1 7 * 1 . o 
0.402 <m> 



Ooe rating Varia:?le« 



InDULt variables are: 
rieat I oac 

Coolant inlet tefiip 
Coolant outlet temp 
Coolant mass flow rate 



50.00 (Ni 
45.00 <C> 

45.72 (C^ 
^.6329E-02 (k'g/s> 



Fluid prooerties evaluated 
Specific neat <Cp> 
Viscosity (hu> 

Thermal cond ( If ) 
^randtl number <Pr> 
Density (Rho) 

Kinematic vis (Nuv) 
Coef ther exp (Beta) 



at 45.36 (O are: 

= 4224.7 (J/kg.K) 

= 5.8661E-04 (N.s/m‘2) 
= 0.6380 (U/m.K) 

= 3.884 

= 990.0 (kg/m'3> 

= 5.9255E-07 (m''2/s> 

= 4.1266E-04 (l/K) 



*** Calculations 



Preliminary calculations: 
Reynolds nunoer (Re) 
t-riction factor (stat) 
Heat flux (Cipd) 

Mean ^luid vel (Vm) 

Rpw 



7 . -^141 E + 03 
3,^02E-02 
4.060E+03 (H/m'2) 
9.26E~01 (m/s) 
3600.00 



Resul ts : 

Nudb 
49.52 
Ka 

7.39E-^03 

RaRe 

5.50E+07 



N unak 
55.30 



^lunf Nuuim 

ni.80 31.56 



KO 

5.16E-0^ 

7.506E-05 



Ra-^Re*^r 
2. 14E + 08 
Tv 

0,026 



Mus Nuwm2 
^4.84 102.34 

•it Nak 1 

1.63E-02 4ZI.90 
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APPENDIX C 

THERSOSIFHCN AHAIYSIS {COHPOTSR PROGRAM AND RESULTS) 



The following program was written to give the results of 
the correlation for the heat transfer for the rotating, 
closed-loop ther mosy phon discussed in Chapter III. It was 
written in HP-Basic for the HP-9826 computer. Sub- programs 

for the graphics and the thermophysical properties are not 
shown. A sample run is included immediately following the 
listing. 

-Pages IIU-117 Program Listing 
-Page 118 Sample Calculation 



1 1 4 



•000! 




THSY 






• 1 . 


t tll.J ■ 


eDur ar V 






* vjy ! 










1030 


COM /Cnin/ Xri i n , Y.ti i n . Sr x . Sf v 






1040 


PRINTER IS 1 








;050 


BEEP 








■ oso 


PRIMT USING • 


2X . Default vaiues: 






1070 


D=4.763E‘3 


t Tube ciameter (m) 






'080 


^=.8230 


1 'ube length ( n ) 






1090 


R=.4G21 


! Radius of rotor (m) 






* ^00 


Kff = Q 


! Friction Factor for 


bencs 




1 ' 10 










1 120 


PRINT USING • 


•4X,-"TuDe diameter 


--.Z.4DE.’- (m)' 




1130 


PRINT USING • 


•4X.-"TuDe lengtn 


— .DD.D/'" (m)-' 




I '40 


^RINT USING 


'4X/*"CI<t Length (Lt> = 


-••.DD.D.-- (m)-' 


-- :Lt 


1150 


PRINT USING 


'4X. --Radius of rotor = 


•-.Z.3D.-- (n)- 




1 160 


BEEP 








1170 


INPUT "OK TO 


ACCEPT DEFAULT VALUES 


( l=Y.0 = N)?’Md 




1180 


BEEP 








i ;90 


INPUT ’’LIKE 


=1 HPRD COPY ( 1=Y.0 = N)7- 


. Ihc 




1200 


IF Ihc=1 THEN PRINTER IS 701 






1210 


IF Id=1 THEN 


1280 






1220 


BEEP 








1230 


INPUT -ENTER 


TUBE DIPMETER (m)".D 






12^0 


BEEP 








1250 


INPUT -ENTER 


TUBE LENGTH <m)".L 






1250 


BEEP 








1270 


INPUT "ENTER 


RADIUS OF ROTOR <m)-.R 






1280 


PRINT USING ' 


"lOX,--**^ Geometric Variables 




1290 


PRINT USING ■ 


’• 1 4X . -"TijiDe diameter (D) = --.Z.3DE.-*’ 


( m ) -' 


1300 


PRINT USING ' 


-14X.--TuOe lengtn (L> 


= -••.DD*DD.-*’ 


in) 


•310 


^RINT USING ‘ 


'• 1 4X . -"Clf t lengtn (Lt) 


= --.DD.DD.-- 


( m ) 


1320 


PRINT USING 


-14X.--L/D 


= •-.5D.D-:L/D 


, 


' 330 


^RINT LfSING ' 


- 1 4X . '•"Rotor radius (R) 


= --.Z.OD.-" (n) : 


^340 


PRINT 








■ 350 


PRINT USING 


- 1 OX . -"ComDuted values 


are : 




1360 


PRINTER IS 1 








1 370 


BEEP 








‘380 


PRINT USING 


"2X. -••SELECT OPTION: 






1390 


PRINT USING 


-4X.--1 Single ooint 






1400 


PRINT USING 


"4X."”2 Na versus RPM" 






1410 


IN^UT Op 








1420 


IF Uo>l THEN 








1430 


BEEP 








14Z|0 


J;*0 








1450 


printer is ' 








•460 


INPUT "LIKE 


TO PLOT ( l = Y.0=N)?'’.0kplot 




1470 


BEEP 








•480 


INPUT "Nadb- 


1 . Nunak 1=2 . Nuum=3 . Nus=4 . Nui =5" . Ian 




1490 


IF OHplot'1 


THEN 






1500 


CALL Plot(Ji.X.Y.Dt.Tc> 






1510 


INPUT -LIKE 


A LABEL1’(1-Y.0=N)".I1 






1520 


IF 11=1 THEN 


CALL Lauel 






1530 


END IF 








1540 


END IF 








1550 


Nstep= 1 1 








^560 


IF 0p=2 THEN 








•570 


BEEP 








’580 


INPUT "ENTER 


RPM RANGE <NIM ..NAX)" .Romi .Rpmn 
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'59G R^>m = Rpfn^ 

0 0 _ P L 

;610 BEEP 

*E2u I»^»ruT ‘’EMTE-^ I'f- -('ll mK‘\ 

1830 END IF 
BEEP 

1650 IF Ihc=1 THEN PRINTER IS 701 
1660 IF lnc=0 THEN PRINTER IS 1 
1670 f^RINT 

1680 PRINT USING “lOX.””-^** Operatina Variables 

^590 ^RINT 

•700 PRINT USING " 1 OX . “’’Input v^ariar>ies are: 

1710 BEEP 

1720 INPUT “ENTER HEAT LOAD <H)“.Q 

1730 PRINT USING “MX.’“’Heat load - . HD . DD . <W) : 0 

17HC BEEP 

1750 INPUT “ENTER COLD SIDE TEMP <C)“.Tc 

1760 PRINT USING “1HX,““Coid side temp = “”.HD.DD,““ <C> :Tc 

*770 Dt=5 
1780 Th=Tc+Dt 
1790 Ta= ( Tc-^Th) * . 5 
1800 Dtr=Dt/L 
1810 i^RINT 
1820 CD*FNCPu(Ta> 

*330 Mu=FNMmj(Ta) 

18H0 K-FNKui(Ta) 

1850 Pr=FNPru<Ta) 

1860 Rho*FNRhouj( Ta) 

1870 NijLv/=MuL/'Rho 
1880 Beta=FNBeta( Ta) 

1390 0mega=i\Dm*2*PI/60 
1900 7h=Tc+Dt 
'910 Ta=Tc-^Dt/2 

^920 printer is * 

•930 PRINT “TA = ”.Ta 
1 9H0 Rhoh-FNRnoui( Th ) 

1350 ^RINT “RHOH*“.Rhon 
1960 Rhoc = FNRhouj( Tc ) 

1970 PRINT “RHOC=“.Rhoc 
1980 Rhom='FNRho»i.i< Ta) 

1990 Drho=Rhoc-Rhon 
2000 IF Drho<0 THEN GOTO 31 HO 
2010 PRINT “dro”“,Drho 
2020 Nu*FNNuu)< Ta) 

2930 Cp=FNCDu(Ta) 

20H0 Be t a=FNBe ta< Ta) 

2050 Pr'FNPruiaa) 

2060 Nuv=Mu/Rhorri 

2070 PRINT “MU. CP =“.Mu.Cp 

2080 Mf=0/(Cp*Dt) 

2090 Re*H/PI*Mf/(Mu*D) 

2100 Gr = R*(0mega^2>*Beta^( <D/2)'‘H>-Dtt/Nuv"2 

2110 Ra=Gr*Pr 

2120 Prod = Ra*Re*-Pr 

2130 Nuum-.262*Prod^ . 173*H8/1 1 

21H0 D twp=0/ (K‘^Nuum*PI*L ) 

2*50 PRINT “DTWP =*’.Dtup 
2160 Dtu)=Dtuip 

2*70 Grr = R*0mega''2-^Be ta*D 3*Dtu/Nuv'2 

2180 Nu=Pr/Grr*L t /R*D/L*RNorri/Rhoh* (KfT + (256/Re^Lt/D) )*Re' 3 
2^90 Dtuc-0/(PI*L-K:-Nu) 
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2200 ABS ( D tu<'~Dtu ^ ^ 0 

22^0 IF Rpn>30C0 THEN 

^woO ou » u . ->oo 

22A0 END IF 

2250 IF Rpm>1000 THEN 

2250 Dt-Dt-^1 

2270 GOTO ^900 

2280 END IF 

2290 IF RDm>500 THEN 

2300 Dt=Dt+1 .5 

23:0 GOTO tSOO 

2320 END IF 

2330 IF RDiTi> = 300 THEN 

2340 Dt=Dt^4 

2350 GOTO ;900 

2360 END IF 

2370 END IF 

2380 Vm=4^Mf/(Rhom^PI-D''2^ 

2390 Ro=Vm/( 0meqa*D > 

2400 Jay=Omega*D *2/Nav !From STEPH-Ro tat lonal Re 
2410 Dtt=<Th-Tc)/L 

2420 Gr = R* (Cmega"2 ) ♦Beta-^ < < D/2 ) 4 ) *Dt t/Nuv '2 

2430 Ra=Gr*Pr 

2^40! 

2450! Laninar Value 
2460 Nulan=48/11 

2470! 

2480! Turbulent Value f ron D 1 1 tus-3oe I ter 

2490 Nudb-. 023-Re .8-Pr' .4 

2500! 

25’0! Stephenson Correlation 

2520 Nus=(Pr . 4/ . 367 ) - . 007 ! -Re .88-jay'.023 ! A i r . TurD-EQN2 . 23-cor recred for uat 

o »• 

2530? 

2540! Nakayama Correlation for larunar conditions 
2550 J2=J-2 

2560 Sq = 2/P-< P$0R( P77/4-( i/Pr ‘2)) > 

2570 Cf = 1 - . 486- < ( 3-Sq~ 1 ) " . 4/ ( Sc- < Sq-Pr-SQR ( 5 ) -2 ) ) ) - J2/ ( Ra-Re ) .6 
2580 .Nunax 1=48/; 1-. 131/$a-(3-Sq-l . 2-( Ra-Re ) ‘.2-!/<* + <t/10-Sq-Pr>> 

2590 Nus = Mu !Temp Value Nus'S'Nijl 
2600 Prod=Ra-Re-Pr 
2610! 

2620! Hoods-Morris Correlation for Lajninar 
2630 Nuu/m= . 262-Prod" . 1 73*48/1 1 ! Eqn 2.22-WM3 
2640 Jl=Jay/8 
2650! 

2660! Woods-Morris Correlation for Radial Pipes 

2670 Nuuim2= . 0 1 5-Re" . 78- J 1 " . 25 Ifrom um2 

2680 IF lhc=0 THEN 

2690 PRINTER IS 1 

2700 ELSE 

2710 PRINTER IS 701 

2720 END IF 

2730 IF Jj=0 THEN 

2740 IF Jj = 0 THEN PRINT USING *M OX /•"Resul ts : 

2750 PRINT USING ”1 OX , ”'*F luid properties evaluated at . DDD . DD , (C) are: : 

T a 

2760 ="RINT USING ’* 1 4X . ’’’’Spec i f ic heat <Cp> = •*’*. 4D . DE . •**’ ( J/k g . K ) :Cp 

2770 PRINT USING ”1 4X . ’‘’’V iscosi ty <Mu) = Z . 4DE . (N.s/m'2) :Mu 
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2780 




USING 


” : 4X . 


’’“^Kernal cond < k 


) 


= 


t t » 1 


. Z . < .*i/fri . k > 


2790 


=^RINT 


USIKG 


”'4>:. 


*”'Prandii nujnoer 


(^T ) 


= 


MM 


. Z . 3DE" : ^'T 


2800 


PRINT 


USING 


”14X. 


“’’Density <Rho> 






»l tt 


. 4D . 3 . “ “ 9 / m 3 


■*' J * 


" 21 *. 


U C I L C 


• • ♦ . vy 

I ^ y. 1 


”’*k i rstrna 1 * 


> 


- 




. Z . HZZ . V F ' 2 '3 


woZG 


- rv iin 


USinb 


’ < . 


therm exD<Deta> 


= 




. i . 4UL . 


2830 


PRINT 


USING 


”14X, 


’’’’Reynolds mijnber 


(Re) 


= 


9 9 9 % 


.Z.3DE“:Re 


2340 


PRINT 


USING 


”14X. 


'”’Nean fluid vel 


(Mra) 


= 


It ft 


.Z.2DE,““ (n/s)“ 


2850 


PRINT 


USING 


”14X, 


, ““Mass floLi rate 




= 


MM 


.Z.2DE.’”’ <kg/s> 


2860 


PRINT 


USING 


”i4X, 


, “’’Temperature (hotside) 


= 


f« If 


.DDD.DD.’”*(C) 


2870 


PRINT 


USING 


”14X. 






* 


99 t« 


, DDDDD . . : Rpm 



: ^"0 
' : Ci9ta 
' : V m 



ilh 



2880 
2890 
2900 
2910 
2920 
2930 
2^^0 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3 *00 
31 10 
3'2C 
3^30 



END IF 
PRINT 

IF 2Wp.oz = 0 CR (Gkpiot = ' "HEN 

PRINT USING *’10X.’*" Niidb Nanair i Nuum 
PRINT USING ••:3X.6(3D.DD.2X)**;Nudb.Nunaki .Nuum.Nus 
PRINT USING ’M^X.”” Ra Ro Ra-Re-Pr 

PRINT USING ’M4X.4(Z.2DE.2X) ,2D.DD”:Ra.Ro.Prod.f1f 
END IF 

IF 0<pLot=1 THEN 
IF Ian=1 THEN Y=Nadb 
IF Ian=2 THEN Y^Nunakl 
IF Ian=3 THEN Y^Nuuin 
IF Ian=4 THEN Y = Nu.s 
IF Ian-5 THEN Y=<48/1 1 ) 

X = Rpm 
J T * J j + 1 

CALL PlotCJj.X.Y.Dt .Tc) 

PRINTER IS i 
END IF 

IF 0p=2 THEM 
Dt = 5 

RDm=RDn-^ 1 0 . 1 " 

IF RpmXRpmh THEM 3130 
GOTO 1890 
END IF 
BEEP 

INPUT ’’ANOTHER RUN < - Y , 0-N ) P” . I r 



Mas 



Mf 



3150 PRINT ”PU” 

3.160 IF THEN 1360 

3170 INPUT ’’WANT TO LABELP( 1 -Y , 0-N ) ” . 1 1 

3180 IF 11 = 1 -THEN CALL ^aoei 

3190 END 
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^ J cr MCf t r 1 C V CL 4 * J j. c* i 

TuLDe diarieter ^i)> = 
length 

Ck t length (Lt > 

L/D 

Rotor radius <R) 
ComDuted values are: 



^ . / b3£~ 0? < o > 
.32 (n) 

2.-^5 Cn) 

172.8 
0 . ^02 ( m ) 



ijperaring Var.aoies 



InoLLt variables are: 
reat load 
Cola siae remp 



50.00 (U) 

^5.00 (C) 



Results: 



- 1 



evaluated at S2 . 75 (C> are: 



uid prooerties 

Spec 1 f 1 C neat ( Cp ) 

1SC03 1 ty ( Mu^ 

Thermal cond < k ) 
Prandti numoer (Pr> 
Density (Rno> 
Kinematic vis (Nuv) 
Coef therm exp<Beta) 
Reynolds numoer ^Re) 
Mean fluid vel (Vm) 
Mass flow rate 
Temperature (horside) 
RPM 



= ^310. 1E+Q0( J/k 
= -^.^1 72E-0^ (N. 

= 0 . 0 ^ 0 ^ (W/m.k > 

= 2.S02E+00 
= 989.1 9kg/n'3> 

= 4.4991E-07 <fi'2/S) 
= 5.39^0E-0^ n/K) 

= I .978E^02 
= 1 .87E-02 (m/s) 

= 3.27E-04 (kq/s) 

= 80.50(0 

= 3600. 



Nudo Nunai' 1 Nuum Nus 

2.42 5 1 .30 34 . 1 ■ .65 

Ra Ro Ra'^Re-^^r 

S.13E^05 :.0^E~02 3.52E^08 3.27E-0^ 



119 



ot 01 



A£P^DII D 

CLOSED TSO-rH&SE THE E20S YPH ON SHALYSIS (CO.^POTBH P50GRAH AND 

HESOLTS) 

The following program was written to give the results of 
the correlations for the sonic, boiling, and condensing 
limits as discussed in Chapter IV. It was written in 
HP-Basic for the HP-9826 computer. Sub-programs for the 
graphics and the t her mophys ical properties are not shown. A 
sample run is included immediately following the listing. 
-Pages 120-124 Program Listing 
-Pages 125-127 Sample Calculation 
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I 





"I,£ 


~ 








: *‘’ar 


- r, F . ' ' j - ; 






. ^ A • 










^030! 


2170.CAu:S; «000. 


^Ll^: 5260. FNS: 61 OC 


. uABEL : 




’ : 








•050 


CC- /Cam/ <n i n . ' 


1 : n . S f X . S f y 






*060 


PRINTER IS • 






’ 070 


BEE? 










Dp^SjT ilST'^K^ "oy " 


"OefajLlt values: 






lOrO 


D-^.763£-:- ! T'loe Oiaaeter <m> 






’ «00 

1 < 


■-€■=.3230/2 I : 


!jJDe I encrri ( Evap M a > 






! • :o 


— N = ' -■'* i"" 








* ^30 


5=,j02’ * - 


ao i L'.s o - r o 1* o r ^ a ^ 






• ’ -20 


=’-INT JSING ”RX/’ 


“T:xos 0 lane re r 


. Z . 4 D - . ( n ) 


• ; [) 


^ *50 


'PINT USING ••-•/,'• 


lengr n < E vao ) = 






' ■ 60 


PRI-r USING "-x.” 


*'^oi:'e Lena th ( Cone ) = 


. DD . J . ^ m > 


Lc 


1 !7Q 


PRINT USING ••-ix/’ 


'‘Paeius of rotor 


■*“.Z.3D.““ (n)' 


A 


■ ^ so 


3E£^ 








1 *30 


INPUT "GK 70 f^iCCEPT DEFAULT UfiLUES < ! 


-y,0=N)l'“.Ic 




^200 


3E£=> 








_ j 


INPUT ••'.IKE P hPRD ISP'-' ( ■ --'.■3=-N>''". I 


he 




* "20 


* r — 1 ^ , 








•230 


PRINTER IS 701 








*2^0 


ELSE 








; 250 


PRINTER IS • 








1260 


END IP 








i27Q 


E^ -G=l then i3b0 








*230 


BEEP 








;230 


INPUT “ENTER TUBE 


DIAN£"ER <a>**.D 






1300 


BEEP 








I 0*0 


INPUT “ENTER TUBE 


LLNG l^H ^ E vap ) (in)'‘.Le 




1320 


BEEf^ 








1330 


IMPU^ “ENTER TUBE 


-ENGTH(Ccna) (n))“.L 


c 




1 3^0 


BEE^ 








I 350 


INPU^ “ENTER RPDI 


US GF RO 1 GR < n > “ , R 






^ 360 


PRINT USING “*0X. 


Geome:ric Van 


adles --- 




!370 


PRINT -USING “l^X. 


’*‘*TaDe aiajneter <D) 


= "“.Z.3DE."” (rr.)' 




1330 


PRINT USING *’l/iX. 


““TiiDe i eng 1 6 ( Evap '> 


(Le> = "••.Z.3D."‘ 


* (rn) 


1390 


=^RINT USING “l^X. 


ieng tn < Cond > 


<Lc) = “”.Z.3D.““ (n> 


1^00 


PRIN^ USING *'i^X. 


’’“Rotor racius CR) 


= .Z.3D."’ 


( m > 


MIO 


PRINT 








1^20 


PRINTER IS 1 








M30 


BEEP 








I^^O 


=*RINT USING “2X.“ 


“SELECT OPTION: 






1^50 


PRINT USING “^X.“ 


1 Q versus /.Fill 






M60 


PRINT USING “4X.“ 


797 ??997?9??'*’*** 






1 ^70 


INPUT Op 








1 ^80 


iF Gp > 1 1 HEN 








1^90 


3EE^ 








1500 


PRINTER IS 1 








J5I0 


PRINT USING “PX.“ 


“NOT INPLENENTED YE'" 






*520 


GOTO 1^/JO 








1530 


ELSE 








15^0 


BEE^ 








1550 


G=0 








*560 


Jt = 0 








1570 


D ^ r = * 








’530 


INPU^ “LIKE TO PL 


GT < 1 = Y . 0=N ) P“ . uk p lo t 





: Lc 
;R 
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Jr>oO 'j ~ ^«J! J , _ = iVSn- , j = ~ JU i * . ~y p«? 

* ^ n ^ M ^ 

:S50 END IF 
IRSO N;:?c='0f' 

IS70 BEE? 

'530 INPUT -EMTE? jPTGN < I =3 . Z=£ . 3 = 3 . •i=C) 'Mo 
1590 F Io=2 THEN 
I ^OG 3EEP 

• ? • 0 INPUT • -UPLIFTS . E-..AoTER” . le 
i720 END I- 

•730 IF :-c=’ 'FEN --IN'EF 13 70! 

1740 IF lhc=0 THEN PRIMTEF IS ’ 

’750 3EE^ 

*7^0 INP’!^ ’*-NT^R ~sa.t <C) ".~sa‘ 

17^0 PRINT USING •*i4X/*“T.5aT<o'’ = . DDD . 0” : Tsa: 

: 730 3EE=^ 

1730 INPy-^ ‘’ENTER RP^“.Rd.i 

^300 =‘RINT USING ” * 4X . 

‘310 BEE’- 

*•320 INPUT "ENTER CGMD. TE;iPERAUiRE" . Ty 

1330 PRINT USING "‘4X.""T ua.1 1 - . DDD . D" : Tw 

;340 At=^I-0 2/4 

1350 *sark=T3at-273. 15 

^ 360 vIp = ^NCDfji "sai: > 

1 870 Nal y ( T-sat ) /? > 

! 380 ^av-^NNav ( Tsat ) 

1890 K=-NKii.'(Tsat ) 

1300 ‘^r--NPrrj( Tsat) 

191Q Rho 1 =FNRhow ( Tsat ) 

*920 H+ g = - iNHf g ( Tsat ) 

1930 Be ra==FNBe tac Tsa: > 

*94Q S iona=»^NS 1 gna< 'sat ) 

1950 Rhov= ' /FNVvs I < Tsar ) 

1 3E0 Nllv = h»j../ /R hov 
1970 Rsn=1.327 

1 980 Cv = CD/Rsn 
^990 Rgc= 4G1.52 ! (J/kg.K) 



2000 


Onega= 


pDn*Z* 


PI/60 












2010 


Acce 1 = 


'Omecia 


2^R 












2020 


IF J; = 


0 THEN 














2030 


PRIM" 
















2040 


^RINT 


USING 


•’;ox. 


""Fluid orooerties evaluated 


at "".ODD 


.DD. "" (C) are;"’ 


I . sar 
2050 


^RINT 


USING 


"*4X. 


""SDecitic dear 


= ’* 


".4D.D."" 


( J/'<g.K)" 


: Cp 


2060 


PRINT 


USING 


"14X . 


""Rario of Scecific Hear 


s= " 


".D.DDD":R 


s n 




2070 


PRINT 


USING 


"!4X. 


""Viscosity (Mu) 


= •' 


•M.4DE."" 


( N . s / m 2 


) ; Nui 


2030 


PRINT 


USING 


"14X. 


""Thermal cone < \c ) 


= ” 


".Z.4D."" 


(H/m.K)"" 




2030 


PRINT 


USING 


"14X. 


""Prandti numoer (Pr) 


= *’ 


".Z.3D":Pr 






2100 


PRINT 


USING 


"^4X. 


""Entnalpy (Hfg) 


- " 


",Z.4DE";Hfa 




2! !0 


PRINT 


USING 


"i 4X . 


""Density (Phoi) 


= ” 


".4D.D."" 


( X g / m 3 ) 


"" ; Rho i 


2120 


PRINT 


USING 


"14X. 


""Density (Rhov; 


= " 


".Z.4D."" 


( Ic q / m 3 > 


: Rhov 


2’30 


PRINT 


USING 


"MX. 


""Coef rner exo (Beta) 




".Z.4DE."" 


( f/K) 


; Be ta 


2^40 


PRINT 


USING 


"MX. 


""Omega 


= ” 


".4D.D."" 


(i/s> : 


Omega 


Z’50 


^RINT 


USING 


"MX. 


""Surface Tension 


= " 


" . Z . 4DE . "" 


< M/m ) 


; S igma 


2^60 


END : = 
















2 '70! 




















PCF = 


•'C '■ 












. 
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qc * 

::co - JO-:' :/j 

22^ 0 Av=( '^/(-_t*(Rhov-^hrii ) > >-< (Phoi/cRhov-Rhol ) • 

'!' ^?Q » 

2230! SCMIC CPLC’JLRriari 
2240 IF Io = i THEfi 
2250 IF J^=G TMEM 
22S0 PRINT 

2270 USING OX . *'*'Son ic Limit Caicui at i or.s : 

2230 PRINT USING ’MOX.””/: F.il Total Meat 

2290 EMD^IF 

2300 = < *'* / ( l t < Rhc'> ~Pho I ) > -Rho I / ^ Rhov-Rho i ) t ) *S0R ( Rs h^Rgc^ T^a-tk '> 

2310 MT = ijci/HT.- 
2320 X--C- 

2330 Y=Qm 

2340 IF Gn< = 0 THEN GC'^0 2400 



2350 


IF OkDlot=" 


HEN 


2360 


CFL. ^lor(o;. 


Jc . Jc , X . . Type . Cx ) 


2370 


ELSE 




2330 


^^I.NT JSIMG •• 


•2:(.3D.3X.2.3DE”:?cP.0m 


2390 


END IF 




2400 


^c->=’00 T 


HEN GOTO 3730 


2410 


IF Ok plot*" T 


HtN 


2420 


Pc‘ - i 0 . 02 


2430 


ELSE 




24^0 


=Pcf - ; 




2450 


J T = , J r + 1 




2460 


END 




2470 


GC'O 2200 




2480 


END IF 




24qo ’ 






2500! 


entrhimment C 


ALCULATICNS 


2510 


IF Io=2 THEN 




^5^0 » 






2530 


IF J;=0 ”HEN 




2540 


PR IN’ 




2550 


PRINT USING •• 


" OX . ""Ent ram ing Limit Clacuiations 


2560 


PRINT USING •' 


lOX. */ r - 1 i 1 Total Heat 


2570 


END IP 




2530’ 


^ALlIS CORRELATION J-=.25 


2590 


IF Ie=" then 




2600 


IF Av=0 THEN 


GOTO 2990 


2610 


Jv=.25/Rhov . 


5* < Omega' 2^R^D*<Rho:-Rhov>) .5 


2620 


Qm = Rhov'* Jv-*A t 


,^HPg 


2630 


END IF 





2640! 

2650! JOiSTZR CGRRELATIGN 
2660 IF Ie=2 THEN 
2670 IF Pc*=0 'HEN ;^cf = l 
2680 Tl = (l/(i-(Pcf/in0)))-^* 

2630! T;-< V(Pcf/100>)-! 

2700 T2=Tl-(Rhol/Rhov)^(2/3) 

2710 IF ^\<0 them 2990 

2720 Xx2=nial/Mav) ' .25-T^ ! . 75^ < Rho 1 /Rhov ) ' ( < 2/3 ) - ! . 75- D 
2730 Xx*V<i+T2) 

2740’ <x=.5 

2750! <x2== < Hul /Huv ) ' . 25*^hov/Rho 1 
2760 IF G=0 THEN G=' 

2"70 Rel-< ’ -Xx)-<G-D/Hai) 

2780 Rev = Xy-G-D/Hiiv 
2^90 “v=.079/Rpv .25 
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2?00 


IF Fe<=‘25C T 


H - N j c = S Q F C F e i / J 


-'O » 

— VJ 


IF Fe>’25-j ''H 


EN Io= . J53 ; 


2320 


Nim = * ( 2 ^ 


* < ^ho 1 -^hcv > " . ' 


2S3G 


Den= ( 1 '»■ < <2 ' . 2 


5> 7^; 


23^0 


Gc = < Niixn/Den > 


C/3) 


2350 


IF nBS<(G-Gci 


/Gc>>.0! ihiiN 


2360 


G = < G ^ ' j c ) / 2 




2370 


GG^G 2760 




2830 


END IF 




2890 


Qrn = V ♦^-jf 


g 


2900 


-ND 




2910 






2920 


Y = Qn 




2930 


IF Qn<-0 THEN 


GOTO 2990 


2940 


IF u1cpiot = ^ T 


HEN 


2950 


C.9LL 


Jc . Jc . X , Y . : ype . ( 


2960 


C! Cp 


— 7 J 


^?INT ..SING •’ 


12X.3D.3X,Z.4DE, 


298G 


END IF 




2990 


IF Pcf>=i00 ’ 


HEN GOTO 3730 


3000 


IF Gi^Diot = ‘' T 


HEN 


3010 


'=>cf = ^c*^i0 ' .0 


2 


3020 


z' cr 




3030 


Fcf = ’^cf 1 0 




3040 


END I- 




2050 


J ’ r - ? 




TOGO 


iF Ie-2 ^HEN 


GGTC 2580 


3070 


IF Ie=i THEM 


GOTO 2200 


3080 


END IF 




3090! 






3!00! 


BOILING CALCULATIONS 


2*^0 


IF Io=3 THEN 




3:20 


IF Jj=0 THEN 




3130 


prim: 




3140 


^9imt using ” 


’ OX , ”‘*3o lima L 


3150 


USING 


Fill 


3160 


END IF 




3:70 


IF ^cf>‘00 'HEM GCT*3 3720 


3180 


Ab='D*L‘?*ACS ( 1 


- ^ c / 5 0 ) 


3190 


Qr»=Ab* . I 3*RHov . 5*H^a^ < Umeqa 


3200 


X = ^c* 




3210 


f = 0n 




3220 


IF Cn<-0 THEM 


1 3230 


3230 


IF Okpiot=l THEN 


3241] 


CALL ^Lot<oj. 


.^c . ,Ja Tyoe , ( 


3250 


ELSE 




32S0 


=^RINT USING 


'12:<.3D.SX.2.4DE 


3270 


END I" 




3230 


:r ^cf>100 THEN GOTO 3730 


3290 


J j=Jj*1 




3300 


IF GkDlot=! then 


3310 


Pc — lO . 


02 


3320 


ELSE 




3330 


Pcf =^cf 1 




3340 


END IF 




3350 


GOTO 2200 




3360 


END IF 




3370» 






3380’ 


CGNOENSATICN 


calculations 


3390 


IF Io = 4 r{_jPfj 





iiTiit Laicaiat ions : 
Total Hear 



QDP’ 
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3 “^00 = ^E’l 

34*0 

3420 vSI'iG “ : OX , “‘'ConGens inc Lidii Caicaiat ions : 

3^30 PRINT USING ’V.0X.‘’“7. Piil 'loial Heat 

34^0 END IF 
3450 Dt = Tsat-Tij 

3^60 1 ^ . 68- < Cp-D t /Hf g > ) 

3470 Tgc= < Rho 1 * ^ RHo I -Rhov ; ^Gmeqa 2-R-Hf gp-^E ' 3/ ( D*Mul -D t > ) .25 
3480 DEG 

3490 PSi-’80-ACS( I - (Pc‘%'50) > 

3500 IF P->i> = 0 AND Pni< = n0 THEN F = -MF!(Phi) 

35i0 IF Phi>I:0 AND ?hi<=^50 THEN F=FMF2(PhL> 

3520 IF Phi>i50 AND Phi<-i80 THEN F = “pNF 3 < Ph i > 

3530 IF =^hi>!S0 THEN PRINT ”ERRGR . PHI> 1 30 ! ! ! ‘’ 

35-0 RAD 

3550 Qfn=Av^D t-P-Tca 
3560 X=Pc‘ 

3570 '^-Qn 

3580 IP ■jrn<=0 THEN GOTO 3700 

3530 I' GlfPior=l THEN 

3600 CALL p Lot < J j , Jc . Jd . X . Y ♦ Type . Cx ) 

36 1 0 ELSE 

3620 PPINT USING ” 1 2X . 3D . 3X . 3D . D” : Pcf . 0;;, 

3630 END IF 

3640 IF DkpLot=' ^HEN 

3650 Pc^=Pc— '0‘ .02 

3660 ELSE 

3670 Pc‘=Pcf-l 

3630 END IF 

3630 

3700 IF Pc*>=Mjn THEN GQTQ 3730 
3710 GOTO 2200 
3720 END IF 
3730 3EFP 

3740 INPU^ ‘•AiNQTnER nUU < ; = Y . 0=M ) T'M r 
3750 IF Ir=^ THEN 1200 
3750 3EEP 

37T0 INPUT '’HANT 'Q ._ABELT' ( I = f . 0=N ) 'M 1 
3780 IP Il-l THEN cal:. Lacei 
3790 END 
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.eonecric Oar 


. ao - 3 








Our-e dianeier 


^ [■ > 


= 


u , 


763E-.-.- <r 


Ti ) 


1 ! n g« 1 eno t n < E 


ap ) 


(L 


e ) 


^ ^ i 


■ n ) 


TuLoe leng:h(C 


ond ^ 


( L 


.c) 


_ -y * 


< n ) 


Rotor radius 


(R) 




0 . 


402 ( n ) 





at(C)= 30 

^ = 360 C 



•“ X 2 oroDert^e'r evai ■^a''-_-G 
SDeci*ic !“iear < Co ) 



at SC . 03 < 1) are : 
= 1 . S < J /k 9 , y . ) 



Ra^io of Soecific -leat 
Viscosity (Mu> 

TSer'^al cone <k> 
Prandtl nuunDer (Pr> 
Enrhaioy (ri*g) 

Density (Rnol) 

Dens i ty < Rhov > 

Coef ther exp (Beta^ 
Onega 

Surface lens ion 



3= ’ . 327 

^!^336E-0^ (N.s/n 2> 
= 0.6638 U / n . K ) 

= 2.313 
= 2. 30576^06 
- 971.8 ( k q/n 3 > 

= . 3 < I< q / n 3 ) 

= 6.4578E-04 (1/K) 

= 377.3 ( ’ /$) 

= 7.704SE-0Z <M/n> 
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Conaensing Lifrut Caicaiations 
‘/.‘'ill Total heat 



1 


66 . 1 


2 


64 , 7 


3 


S3 . 4 




62.2 


5 


61.1 


6 


60.- 


7 


58. 1 


8 


56.8 


9 


55.5 


10 


54.4 


1 1 


53.2 


12 


52.1 


i 3 


51 .0 


1 A 


45.9 


15 


48.9 


16 


47.9 


17 


46 .9 


18 


45.9 


19 


45.0 


20 


44 . 0 


21 


43.1 


22 


42.2 


23 


4 1 . 3 


2^ 


4C‘.5 


25 


?9 . 6 


26 


38.8 


27 


37.9 


23 


37.' 


"'9 


36,? 


30 


35.5 


31 


?4 . 7 


32 


OO Q 


oo 


33^4 


34 


r 

Ot. . D 


35 


3^ .7 


36 


30.8 


37 


30.0 


38 


29 . 2 


39 


23 .4 


40 


27,6 


41 


26.9 


42 


26,1 


43 


25.4 


44 


24.6 


45 


23.9 


46 


93 2 


47 


22.5 


48 


21 .8 
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43 


*1 ^ 1 


50 


20 *.^ 


'D I 


I 3 .8 


52 


19.1 


53 


^ 8.5 


54 


1 :' 


55 


17.2 


56 


16.6 


57 


! 6.0 


58 


15.4 


59 


14.9 


60 


*‘ 4.3 


6 1 


! 3.7 


62 


13.2 


S 3 


12.6 


64 


12.1 


55 


11.6 


66 


11.1 


67 


10.6 


58 


10.1 


53 


9 . b 


70 


9.1 


71 


8 . 6 


72 


8.2 


73 


7.7 
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